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NOTATION 

X e A means that x is a member of the set A. 

A <^8 means that A implies B. 

{x|x has property a} denotes the set of all x such that x has 

property A. 

G:A b means that the operator (or function) G maps the set A into 
the set B. 

V denotes . 

* dt 

Re z denotes the real part of the complex number z. 

Im z denotes the imaginary part of z. 

X < * means that x is finite. 

x(t) = a means that.x(t) = a for all t. 


Lim x(t) * a , or x(t) -*• a as t 

t-*® 


means that for all n > 0 there 
is a T such that |x(t)-a|<n 
for all t ^ T. 


sup X denotes the supremum (or least upper bound) of the set of 

xcA numbers A, i.e. the least number y such that x < y for all x C A. 

inf X denotes the infimum (or greatest lower bound) of the set of 

XCA numbers A, i.e. the greatest number z such that x ^ z for all 

X C A. 

stp C is the function on the real line defined by 
(-a,b) j-a, c < 0 

stp a =1 0, cr » 0 

(-a,b) I b, CT > 0 

Sod CT is the function defined by 

(-a»b) Sod o = CT stp 0 . 

(>a,b) (-a,b) 

IR denotes the set of all real munbers. 
j^mxn denotes the set of all mxn real matrices. 

1R^ denotes the set of all n-dimensional real vectors. 


A matrix is denoted A, a column- or row-vector is denoted b. 
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A scalar is a one-dimensional vector, i.e. a real number. 

£ denotes the zero matrix. 

^ denotes the unit matrix. 

A* denotes the transpose of the matrix A. 

tr A denotes the trace of A. 

det A denotes the determinant of A. 

S** ^ denotes the set {x e R^jjc'x * l}. 

A > means that the matrix (A-E is positive definite, i.e. 

~ X* (A-b)x > 0 for all x e Ir", x 0. 

QED denotes the end of a proof. 
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CHAPTER 1 


SWITCHING VOLTAGE REGULATORS 

§1.1 Introduction: Design Objectives 

This chapter addresses Itself to the question of stability in DC to 
DC convertors. In order to better understand the role of the stability 
questioni we shall first review the convertor design task as a whole. 

The designer of a DC to DC convertor usually has in mind four primary 
design objectives, which are (i) Efficiency (ii) Stability (ill) Regula- 
tion (iv) Smoothing. There are secondary considerations such as size, 
weight, and cost, but these can be considered for any one design only when 
the primary objectives have been met. 

In order to meet high efficiency requirements any conversion process 
considered should be inherently lossless, i.e. given ideal components the 
process considered would yield a 100% efficient convertor. We shall there- 
fore restrict ourselves at first to the study- of lossless circuits , auid 
concentrate our efforts on the theoretical problems involved in meeting 
simultaneously the requirements of stability, good regulation, amd good 
smoothing. Ovr components available are inductors, capacitors, transfon;:ers , 
and switches (i.e. ideal transistors, thyristors, and diodes). The only 
resistors involved are the load resistance and possibly the source resis- 
tance. The intentional introduction of resistors in the conversion process 
is disallowed. In practice, of course, the non- ideal nature of the com- 
ponents used will involve small losses. Our justification for ignoring 
these in our initial analyses is that such losses tend to aid in stabiliza- 
tion rather than hinder, so that the most difficult design case occurs when 


- 1 - 



^ 1 
i 

i 


1 


- 2 - 


these losses vanish. In fact the idea that dissipation aids in stabiliza- 
tion is expressed formally in the Positive Operator Theorem^ which wc use 
in Chapter 1 and which we shall discuss further in Qiapter 2. We therefore 
phrase our theoretical problem as "Given an input DC voltage source, 
produce a stable power convertor to deliver power at some other DC voltage 
to a possibly time-varying load, using only inductors, capacitors, 
transformers, and switches, and meeting some prespecified requirements 
on regulation and smoothing." 

There are two types of regulation to consider: 

(i) "Line regulation": to make the output voltage insensitive 

to input voltage variation, 

(ii) "Load regulation": to make the output voltage insensitive 

to load resistance variation. 

As well as being concerned with the steady-state regulating capabilities 
of the convertor, the designer may be interested in ensuring that the 
convertor also has a good transient response, e.g. a quick recovery from 
a temporary short circuit. For a periodically varying load, the output 
impedance of the convertor as a function of frequency may be important. 

Since all DC to DC convertors necessarily involve some AC process, 
([23], [38], [39]), it follows from the smoothing requirement ti^at in 
general a low pass filter must be utilized on the output side of the 
convertor. However a low pass filter has an inherent lag associated 
with it, so that by lowering the cut-off frequency of the output filter 
the d^'isigned may very well degrade the transient response performance. 

For a given low pass filter the transient response can be improved by 
Increasing the closed loop feedback gain, but stability considerations 
impose limits on this gain. 
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If the smoothing requirement is dropped completely, then the problem 
may be easily solvable in the ideal case. For instance, one may have a 
pulse-width-modulated system with the design requirement that the Integral 
of the output voltage over each modulation cycle be constant. Designing 
a modulator to meet this requirement is then easy, in principle. 

To summarize, we see that the basic theoretical problem confronting 
the designer of a DC to DC convertor is that of meeting the three simultan- 
eous requirements of stability, regulation, and smoothing, while using 
only lossless components. Further, these three requirements cannot be 
considered independently, and in order to obtain an optimiim compromise 
between the demands of good regulation and good smoothing, a designer must 
thoroughly understand the stability problem. 

§1.2 A Classification of Convertor Types 

The stability question can be tackled only after an extensive 
classification of the different types of convertor, since different 
mathematical methods will be applicable in each case. The first character- 
istic of a convertor is whether its output voltage is to be larger or 
smaller than the input voltage, and the second characteristic is whether 
or not transformers are involved. In this chapter we consider transformer- 
less downconvertors , or switching voltage regulators. Even in this class 
there are many different analysis situations, depending for example on 
whether or not a load is present, whether or not the load is purely 
resistive, whether or not the smoothing filter is lossless, whether or 
not the voltage source hais an internal Impedance, whether or not the source 
impedance is purely resistive, and whether or not the source voltage is 
time- varying. 
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5l»3 Switching Voltage Regulators without Source Impedance 

Given a fixed DC source voltage E, we wish to obtain a near~constant 
voltage of aE with 0 < O < 1. The scheme we consider that i.own in 
Fig. 1.1. The low'pass filter is of the form shown in Fig. 1.2. The 
switch represents a suitable interconr.oction of diodes and transistors or 
thyristors. The scalar u is the control variable, which takes on the 
values 0 or 1 depending on the position of the control switch, h load 
resistance may or may not be present at the output of the low pass filter. 



Pig. 1.1 



Pig. 1.2 
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Th« approach we follow here ia a state- feedback approach « that is, 
we assume knowledge of the capacitor voltages and inductor currents, and 
use this to determine the desired position for the switch at any one time. 
We must choose a feedback control law which brings the output voltage 
as close as possible to the desired value OE, and gives the overall system 
good regulatory and smoothing characteristics. The assumption that we 
have full knowledge of the state is not unreasonable, since a fully satis- 
factory performance will not be obtained without it, and in practice it 
should not be too difficult to obtain a good estinate of the state from 
measurements (of the capacitor voltages , for instance) . The state-feedback 
approach is the most reasonable in a situation such as this. To adopt some 
other scheme such as pulse-width-modulation is to arbitrarily ''mpose con- 
straints which can only hinder the attempt to obtain an op^ ' m overall 
performance. 

In steady-ste.te operation it is fairly clear that o switch will be 
working in a periodic way, ruch that the average value of voltage on the 
outi^ut side of the switch is dE. The amount of ripple present at the 
filter output will be determined by the frequency of the switch operation, 
and by the size and number of filter components, since the filter is a 
low-pars one, a higher switching frequency will yield a smaller output 
ripple. However because the particular thyristor or transistor switch 
used will have a finite switching time during which it dissipates some 
power, the operating efficiency decreases as the switching frequency 
increases. We assume that a lower limit on the allowable efficiency is 
prescribed, and thus that an upper limit on the switching frequency is 
obtained. The feedback controller should insure that the switching 
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frequency does not exceeu this allowable limit. Once this limit is pre- 
scribed, the ripple will be determined by the number and size of the 
filter inductors and capacitors. Because of physical bulk it is impor- 
tant to reduce the tc^.al inductance and capacitance involved. For a given 
total inductance and total capacitance, better high-frequency rejection 
for a filter of the type shown in Fig. 1.2 is obtained as the number of 
filter components n incrr;ases. Consider for example the filters shown 
in Figs. 1.3 and 1.4, which have tho same tc'tal inductance and total 
capacitance. 
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Fig. 1.3 
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For Fig. 1.3 

V (.) 

"2<»> - 1 ♦ ^“=‘ 


and for Fig. 1.4 


i ' ' 7 5 3 4 

“ 717) “ ^ ^ ^ ^ ® 


So for a sinusoidal input signal of angular frequency (u 


If 



1, 


H^Cjw) 

hJHw) 


1 - 3LCgj^ + L^C^Ol^ 
1 - 4LC0)^ 



where u 


LC 


H4(jW) 

HjCjU) 
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i.e. the higher-order filter is much more efficient. Thus we are partic- 
ularly interested ;.n high-order filters. Unfortunately, the stability 
question becomes ircreasingly difficult to analyze as n increases. In the 
limit, one might conclude that a transmission line would be the best type 
of filter to use, but difficulties with the stability analysis and with 
the state estimation probably rule out this possibility. 


§1.4 Second-Order Lossless Regulator 
(a) Choice of a F eedback L aw 

We now consider the first member of the series of regulators 
of Figs. 1.1 and 1.2. For the particular control law which we choose 
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the stability analysis is fairly simple. We devote so^ie time to it# 
however# since it turns out that this is the only meober of the series 
for which the analysis is simple# yet it illustrates most of the features 
of the higher-order regulators. Pig. 1.5 shews the system under consider- 
ation. The load resistance R is assumed constant; R « denotes the 
no-load condition. . 
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Letting y, « i/l and y_ ■ V we obtain 
1 2 0 
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We now Introduce time- and amplitude-scaling which allows us to assume 

without loss of generality that E •« L ■» C ■ 1 . For# letting T “ u^t 

and z. ■ — r— we obtain 
1 b 




dT 

9 

ds. 


2 


dt 



-1 



u 


x.e. 


. ~ z “ A z + b u 
dT — — — — 


In what follows we assume rhat ® i» i.e. , t = T. Zj^ and z^ are called 
the state variables emd the vector ^ is called the state vector . 

We want to ntake u a function of z^, (i.e. the control u a function 
of the measured variables z^^ amd z^) so that z^ comes as close as possible 
to its desired value a (where 0 < a < 1 ) . That is , we want 

lim z-(t) e a , 

t -► 00 

for any initial condition (Zj^(O), Z2(0)). We call this global asymptotic 
stability about a . Now if z^ < a one would expect that the switch 

should be in the u = 1 position, and if z^ > tx that the switch should be 
in the u = 0 position. This leads us to try the feedback control 


u 


1 , Z2 < a 


0 , Z 2 > a 


- 4 — 
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(We postpone till part (b) a discussion of what happens at ■ 0.) 

Using the methods of §1.4 (b) , (c) , (d) , or (e) we can show that this 

feedback control does give the desired asymptotic stability about z^ ^ 0( 

provided that R has a finite value. When the load is removed (R ® “) 

asymptotic stability is not obtained. Indeed in §1.9 we give consideration 

to the effect of an inevitable small lag in the feedback controller, and 

find that in this situation a small lag yields instability about z^ = ci; 

in fact an oscillation whose frequency is of the order of ■ is 

2Tr*^ 

obtained. Note that the nO“load case is the most difficult to stabilize, 
and for this reason we shall henceforth analyze only the no-load case. 

The reason that this initial choice of control law is unsatisfactory 
is that it takes no account of the inductor current z^^. In conventional 
servomechanism theory terms, we need to provide some "rate feedback" as 
well as "output feedback". If is slightly less than a while is 
large and positive then it should be fairly clear that we want u « 0 rather 
than u * 1. Since z^ = z^ when R = “ let us consider the control law 


1 , Z 2 + 6*2 ^ 


0 , Z2 + 3z2 > OL 


for some 3 > 0 . 


In §1.4 (b) , (d) and (e) we shall show that this control law gives a 
regulator with good stability about Z2 » ot for all values of R. We discxiss 
implementation of this control in §1.9. The parameter 3 can be chosen by 
the designer to give a good transient response to changes i.n operating 
conditions, such as changes in a, in E, or in R. One value of 3 cannot 
give an optimum response for every condition, so a compromise value will 
have to be chosen. An analog computer simulation shows that setting 3 


between 0.8 and 1.0 gives a reasonable overall performance. 
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(b) Stability by Phase Plane Analysis 

First we make a change of variable so that the desired stability 


is about 0. Let 


Xj . 2Xj 


Xj - 2Zj - 20 


SO that our system is 


» -Xj + 2(u - a) 


Xj-Xi 


the desired equilibrium point being the origin in the (x^,x^) plane. 
Next we introduce the "step" function 


stp x » 
(-a,b) 


-a , X < 0 


X “ 0 


X > 0 


Stp X 

(-a.b) b 


- 12 - 

Bxj + Xj < 0 

Bxj + Xj ■ 0 

Bxj + Xj > 0 
Stp (BX2+X3) 

(-2+20,201) I 


- stp (Bxj+x^) 

(- 2 + 20 , 20 ) 

value of o is Icnown, and abbreviate 

stp to stp. 

(- 2 + 20 , 20 ) 

Some remarks are necessary concerning the use which we shall 
be making of the stp function. This is a discontinuous function, and a 
statement such as 

X « sgn x 

requires rigorous mathematical treatment in order to make strict sense 
for all values of x. However, we shall use the stp function merely for 
its notational and conceptual convenience, since it is the limit of 
various classes of continuous functions such as sat kx (Fig. 1.7) or 
tanh kx (Fig. 1.8) where k is an arbitrarily large positive real number; 
such functions are more accurate descriptions of a “real-world" switching 
function than stp x is. In order to put our treatise on a fim mathemati- 
cal footing we need only replace stp x by sat kx (or isome other suitably 
smooth function) each time it is used, and show that the relevant 


sW 


stp 

(-1,1) 


He now let 


1 , 


u • < a 


i«e. 


“4 


2a - 


so that our system is 


Xi « -Xj 


*2 “ *1 


We assume henceforth that the 
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sot kx 
(-o,b) 

-a. X < 

a b 

kx , " ^ X ^ . 

k k 

b, X >1 
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Fig. 1.8 

conclusions still hold true. In fact, many of our results hold true if 
ctp X is replaced by f(x) where f is any monotone nonlinearity with 



f(0) ® 0. (A monotone function f(0) is one for which f(Xj^) < fCx^) 
whenever Xj^ $ X2*) 

Kow we continue with the phase-plane analysis of 


X + X ® 


-Stp(Px + x) 


Tlie trajectories in the (x,x) plane arc arcs of circles, representing 
piecowisc-simple-harmonic motion. If x 4- Bx > 0 these arcs arc ccntcr(‘d 
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on (O«'-20(), and if x 0x < 0 the arcs are centered on (0,2-2a). Fig. 1.9 
shows the trajectories obtained v^en ^ 0 * 1. The line x *>■ 0x ** 0 

is called the switching line * and we need to consider carefully what happens 
at this line. We have said that the feedback nonlinearity is in practice 
a continuous function, and one way to determine what happens along thr 
switching line would be to replace stp 0 by a function like sat ka, with 
k very large. However beca\ise of the chattering behavior along switching 
lines which is observed in practice, it is better to consider the feedback 
nonlinearity either as a stp function preceded by a small time delay, or 
as a stp function modified to include a small amount of hysteresis. A 
small delay in the feedback path will alvrays be unavoidable in practice, 
representing for instance the switching-time of the power transistor used 
as the control switch. In addition, the designer will want to include 
either a small fixed delay or a small amount of hysteresis, in order to 
limit the switching frequency of the power tremsistor because of 
efficiency considerations, as discussed in §1.1 and §1.3. In this thesis 
we adopt the policy that hysteresis is undesirable, and that a predetermined 
small fixed delay is used. The disadvantage with the hysteresis method 
of frequency limiting is that the switching frequency depends on the load 
resistance, and furthermore this dependence is difficult to emalyze. By 
using the fixed delay method the designer can easily set the switching 
frequency to some desii-ed value which will be independent of the load 
resistance. However, whether a fixed delay or hysteresis or a combination 
of both are employed, the conclusions which we shall reach concerning 
behavior along the switching line are essentially the same. 
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Fig. 1.10 

Typical path for ^c+x ■ -sgn(x+x) 
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We now consider in some detail the behavior of feedback systems 
which have a stp function in the feedback path. For a rigorous treatment 
the reader Is referred to the work of Filippov [13]. Other useful 
references are [1] (Chapters 6 and 12) f [18] (Chapter 6), and [22]. 

Consider the feedback system shown in Fig. 1.11, described by 
the equations 

! i(t) - f (x(t) , u(t) , t) 

u(t) • -stp c x(t) 

(-a,b)“ “ 



(which is a line when n *• 2, a plane when n « 3, etc.) divides the state 

space into two regions, in one of which y > 0, u » -b, and in the other 

y < 0, u a. We shall use the symbol S to denote the switching surface. 

The scalar quantity p * £ x is the distance of the point x from S, with 

• • 

due attention to algebraic sign. The scalar P * £ gives the rate at 

which the point x on a trajectory is approaching S. If we consider an 
arbitrary point on S together with values of p and p nearby we can 
identify several situations with respect to the possible signs of p and p. 
Fig. 1.12 shows the four main cases of interest (we are not considering 
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in practice^ even though the function labelled stp is actually a continuous 
function with bounded slope at the origins if this slope at the origin 
were small enough the chattering would cease. The system is now described 
by / 

I x(t) « f.(x(t) , u(t - T) , t) 

I u(t) “ “ Stp £ :t(tl 

Wo adopt the notation to denote (x, -b, O . i.e. ,£ when £ x_ > 0. 

Similarly denotes ^()t, a« t) , i.e. iu ^ when £ x < 0. We obtain first 
a condition for a state on the switching surface £ x ■> 0 to be an end- 
point. From Pig. 1.12(c) we see that p > 0 implies p < 0 in the vicinity 

* «• 

of i.e. ££ <0. Similarly p < 0 implies p > 0, i.e. £ jf >0. 

Taken together these two conditions arc sufficient for ^ to be an endpoint, 
i.e. £ ^ ^ 0 and £ £ >0. Now we determine the rate at which the state 

chatters along S. Suppose that for t < tj^ the state lies in the region 

£ X < 0, and that x(tj^) ■ (where £ Xq ■ 0) , i.e. the trajectory meets 

S at time tj^. Let Xj * x(tj^ + T) which is the state when u switches from 

a to -b, and let x^^ ■ + be the state when the trajectory next 

intersects S. In the interval [t^, t^ t t] we have u a, so that 

Xj ■ Xq + T f“ + OlT^l 

For t > tj^ + T we have u •* -b, so 

X j ■ Xj + (At - T) + 0((At - T)^) 

- Xq + T(f" - X^) + 0[t^] + OlvAt - T)^]. 

Multiplying this from the left by c and remembering that c x, ■ c x. ■ 0 
wc obtain 
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At - -T 


£(£“-£*') 


c £ 


ignoring the higher order terms in T and At. Thus 


—3 -0 -+ — — 

■ f -(£-£) 


At 




c £ £ - c £ f 


c(£"-f”^) 


Now i£ we were to consider a switching ransition going in the other 
direction we would obtain this same result > even if the switching delay 
was different from T: this can be seen easily by observing that the 

result is unchanged when ^ is exchanged with £ and that the result is 
independent of T. Thus, letting T 0 we see that the system trajectory 
approaches arbitrarily close to the trajectory defined by 

. + - + 

X - f - (f -£ ) 7- 

c (f >f ) 


c £“f^ - £ 

m I ■ ' ■ ■■ I 

c (£■-£■*■) 

That this trajectory remains on £ x « 0 can easily be checked by evaluating 
£ X, %diich is 0. This formula has a simple graphical interpretation. In 
Fig. 1.13 the vectors AB > ^ and AD ■ ^ are the velocities on each side 
of S, at the point A on S. Note that the arrows denoting £^ and f in 
Fig. 1.13 are drawn on the opposite side of 5 to those of Fig. 1.12. The 
resulting velocity vector is x « AC where C is the intersection of BD 
with S. vnien the situation of Fig. 1.13 is as in Fig. 1.14 the resulting 


4 














X ■ £(x»u,t) « A X + b tt 


y « c X 


u " -stp y 
(-a,b) 


£ ■ A X - b b 


f ■ A X + a b 


The condition for a point x on £ ]k ■ 0 to be an endpoint is 


i* e* ^ 


cAx-bcb<0 and cAx-i-acb>0 


-acb<cAx<bcb , 


and the chattering motion is described by 


(£ A £ + a £ b) (A X - b b) “• (£ A x - b £ b) (A X + a b) 

(aH-b)c b 


(c A x)b 


A X - 




A plausible but nonrigorous argument suggesting this result is as follows. 

A trajectoiy on £ £ = 0 satisfies £ £ ® 0, i.e. £ A x - £ b stp £ £ “ 0, 

£ A £ , 

80 Stp £ X can be replaced by Then the equation £ ® A £ “ b stp £ £ 


becomes 




1 


j 


i 











This will define an asymptotically stable motion along S if and only if the 
/ --\ 

matrix £ * IJ. ~ eigenvalues with negative real parts; the 

remaining eigenvalue must be zero since £ F ■ It is an 

interesting and useful fact that the nonzero eigenvalues of F^ are the zeros 
of the numerator polynotdal of G(s) « c(I s - as we now prove: 


Theorem 1.1 


6(8) « c(l s - A) 


^(s) 

p(s) 

n-1 ^ n-2 . 

Vl ^ *•••'^'^0 

.n ^ n-1 . _ n-2 

® ■'■Pn-1® +*“*Pc 


where (A»b,c) is a minimal realization of G(s) and assure that 

p ■ II - *Ph«n e ~ PI a 5SL1SL. 

- \- 


. Then det(I s - F) « 


Proof 


- To 


0 1^0 ‘^r-'Vi^ 


-Pq “Pi 




•nd let 


Then 


so that 



Now (A,b»£) is minimal f and is also minimal, being the standard 

controllable realization of G(s) , ([7 ], section 17), so there exists a 


a matrix P such Lhat 




det(^ E - P) *s det^j^ s - A + ^ 

/ -1 ^ -l\ 

» det(I s - P A.P + P A,P 

-“I- %bi “-1- I 


det P(I s - H) P 



1 



ill if^'f I iti I n 




“1 - 


• det (i ® * H) 

-s-alsl BED. 

Vi 

Theorem 1.1 elso follows from equation (13) of reference (^ !• 

Corollairy i The feedback system of Fig. 1.16 is asymptotically stable in 
the chatteri*>g mode if and only if the numerator polynomial of G(s) » 
£(^8-A) ^ b is strictly Hunritz, (i.e. has its zeros in Re s < 0) . 

A plausible argument suggesting this corollary is as follows: If the stp 

function feedback operator of Fig. 1.16 is replaced by the sat kx function 
of Fig. 1.7, asymptotic stability in a neighborhood of the origin is ob> 
tained if and only if q(s) + ~ p(s) is strictly Huzwitz. Letting k tend 
to infinity makes sat kx approach stp x, and the zeros of q(s) -t- ~ p(s) 
approach the zeros of q(s) , and the result follows. 

Mow vre return to the lossless second order regulator described 


X + X «* “Stp (@x-^x) 

(-2+20, 2a) 


X «* A X - b stp £ X 
~ (-a,b) 


0 -1 j b - j 1 I c - IB 1], X - r X, 


1 0 


a ■ 2 - 20, b » 2o . 


The switching line S is x^ + gXj^ ■ 0, and x is an endpoint if 


-acb<cAx<bcb 


1 


J 


i 



1 


T 


1 


'1 
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whlch becomes 


-3(2-20) < -3xj * 2oB 

or 


-2 + 20 < - Xj < 2o 



Motion along S is governed by x ■ £ x where 

b c 


.-(.-If). 


i.e. 



We note that det(Is-F) » |*(Bs+l ) , in accord with the theorem. Thus 

— — p 

chattering occurs along the switching line x. 4 $x 


-2o/-|— ] < X, < (2-20) 


2 0 in the region 

with notion in this region being detemined 


< Xj < 

\ “ ”(f) *1* 

Three conclusions of practical significance arising from these 


results arei 

(i) That on the switching surface a small value of 6 
is desirable for quick settli-ng (i.e* a short 
transient response ) , though this is not necessarily 
true for the overall transient response, 

(ii) The state vector x will not reach the origin in 
finite time, 


i 


I 





1 


(iii) The tum^on and turn-off delays of the switch 
need not be equal for the foregoing analysis 
to apply. 

An examination of the trajectories of Fig. 1.9 shows that indeed 
the desired stability about (0,0) is obtained, for any initial condition. 

Fig. 1.10 shows a typical path. 

We can see why asymptotic stability is not obtained when 3=0 
by considering Fig. 1.16, in which we see that the switching line x 3x = 0 
is now the x axis. The paths form closed trajectories representing sustained 


oscillations. 



Fig. 1.16 


(c) Stability by Total Gain Linearization 


The system we are considering is 


X + X •* - stp(3x+x) 


I 
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By taking Laplace transfoms we can represent this In the iisual feedback 
system form of Fig. 1.17. Consider now a feedback system with transfer 
function G(s) In the forward path, and nonlinearity f{CJ) In the feedback 



Fig. 1.17 


path, as In Fig. 1.18. Assume that f(0) * 0. We can associate with this 
the linear system of Fig. 1.19, with feedback gain k. 
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Aiserman's Conjecture ([34]/ Chapter 7) states that for the system of 
Fig. 1.18 we would expect to have global asymptotic stability if the 
associated system of Fig. 1.19 is globally asymptotically stable for all 
values of k lying in the range of values taken on by as CT varies 
along the real line. This method of investigating the stability of a 
feedback system is palled the method of total gain linearization . Aizerman's 
Conjecture is not true in general, though it appe«u:s to be true for the 
type of systems considered in this thesis, and it cem be shown to be true 
for second-order systems ( [34] Chapter 7) . The nonlinearity stp O lies in 
the first and third quadrants, with taking on all values between 0 

and 4«. The total gain linearization stability conditions are thus fulfilled 
if the Nyquist locus of G(s) does not rross the negative real axis. This 
is so for the system of Fig. 1.17 for which the Nyquist locus is as shown 
in Fig. 1.20. 

The Circle Criterion is not useful here, oecause the interior 
of the disc is the left half plane, and the Nyquist locus enters this 
region. The Popov Criterion does prove stability for this second order 
system, however we shall not consider its application because it is a 
special case of the Positive Operator Theorem method we consider next: 
the Popov Criterion makes use of a first-order multiplier, which is useful 
for our purposes only for the second-order regulator. 


(d) Stability by the Positive Operator Theorem 

We now make use of some concepts from the theory of Positive 
(or Dissipative) Operators. The most important result which we use from 
this theory is called the Positive Operator Theorem. It is felt that one 
of the main contributions of this thesis is in showing the utility of the 




I 
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ideas associated with this theorem. We devote more attention to positive 

operators in Chapter 2. For our purposes here* a positive operator is an 

operator with input u(t) and output y(t ) , where 0 i t < for which 
/T 

I uy dt ^ 0 for all T ^ 0, and for which y(t) ■ 0 for all t whenever 
^0 

u(t) "0 for all t. This last requirement can be written GO * 0. If the 
operator is a convolutior operator, and if it can be represented by a 
rational transfer function G(s) , then it can be shown (P6] Theorem 1,[41], 
[21]) that positivity is equivalent to the requirement that G(s) have no 
poles in the right half plane, that any poles on the imaginary axis be 
simple with real positive residues, and that Re G(ju) ^ 0 for all co, i.e. 
the Nyquist locus of G(s) must lie entirely in the right half plane. An 
equivalent reqvdrement ([ 2 ], [16]) is that Re G(jw) ^ 0 and if G(s) *» 
then p(s) + q(s) must be strictly Hurwitz, (i.e. all its zeros must lie 
in Re s < 0) . Such functions are called positive real and play an impor- 
tant role in electrical network synthesis, since the driving point Impedance 
of a linear passive network is positive real, euid emy positive real function 
is the driving point impedance of an RLC network. 

nie Positive Operator Theorem gives a sufficient condition for 
input-output stability of the feedback system of Fig. 1.21 in which 
'’2'^2''^1''^2 functions of t for 0 ^ t < *. 



Fig. 1.21 






I I 
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The functions and V 2 are the inputs > and can be used to represent driving 
functions, driving noise, or initial condition responses. By choosing v^ 
and appropriately we can use the Positive Operator Theorem to obtain 
conclusions about the behavior of u, ,y, ,u»,y^ as t To do this we 

need a measure of the "size" of a function of time x(t) } we use the L^^norm 


||x(t)|| 


/• \l 


Using this norm we can only handle functions x(t) for which I x^ dt is 
finite; the set of all such functions is denoted L 2 « There are many func- 
tions of interest «hich are not in L 2 , for example the constant functions. 
This difficulty is overcome by using truncated functions, that is, func- 
tions which are zero after some time T. Now it can be shown thiit for any 
function x(i.) e L 2 for which x(t) is bounded or square-integrable (l.e. 


f: 


X dt < ») #we must have Lim x(t) « 0. An operator G is said to be 

“ lloll 


bounded if there exists a real number M such that 


|x| I 


< M for all X. 


Ai, operator G with input u(t) and output y(t) is said to be strictly posi- 
tive if I uy dt ^ nl u^ dt for some Ti> o, and GO ■ 0. (Alternatively 

Jq h 

we can say that G is strictly positive if G - n X is positive for some 
n> 0, I being the identity operator). 

We can now state the Positive Operator Theorem ((42], [34] 

Chapter 4) t 


Theorem 1.2 

If G^ and G 2 are positive with one of them being strictly positive and 
bounded, then ^2 ^^^'^^ver v^ and V 2 are in L 2 , and 

there exist positive constants and P 2 such that 


i 
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Before applying the theorem we need to consider the requirement 
GO ■ 0 when G is a convolution operator. A convolution operator G(s} 
mapping u into y as depicted in Fig. 1.22 is usually taken to mean the 
operator with input u(t) and output y(t) , 0 < t < *", defined by the dif- 
ferential equation 

p(D) y(t) ■ q(D) u(t) %diere P “ ^ 

i.e. 

Pn y^”^ Pq “ Sn +...♦ q^ u(t) , 

together with a given set of initial conditions 

y(0) , y^^No) ,..., y^"“^No) - 

U(t) 


Pig. 1.22 



If m ^ n this means that we can describe G in state space form by 


x(t) ■ A x(t) + b u(t) 


■ c x(t) + d u(t) 


x(0) - Xq 


where x(t) € IR**, ^ ^ “ G(s ) , and (A,^,£,d) is assumed to be a 
minimal realization. Now the operator G defined thus has the property 
GO ■ 0 if and only if « £. For, G(u(t)) is given t ^ I by 



y"” T ^ ^ 1 — ' T ' " " .j- 

I 1 , , I . ' 
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y(t) - £ Xq ♦ 




A(t-T) 


b u(T) dr + d u(t) , 


thus 


tt(t) H 0 


y(t) - c Xq 


^ - 0 if y(t) H 


Further, 


X * 0 
=0 - 


ix. 

dt^ 


t»0 


Since (A,c) is observable 
since (A,b,c,d) is tainietal. 


for 0 < k 4 n-1, since (A,£) is observable. 

In what follows we shall depict the initial conditions associ- 
ated with G explicitly by means of an arrow, as in Fig. 1.23. If no such 


U 



6 ($) 


Fig. 1.23 


y 


arrow is shown we shall mean that the initial conditions are unspecified. 

Now an operator G(s) with initial conditions as in Fi .1.23 

can be represented as an operator with zero initial conditions followed by 

the addition of an external signal which is the initial condition response 

of G to as depicted in Fig. 1.24. We shall make u<.e of this equivalence 

in applying the Positive Operator Theorem. 

Before applying the Theorem to our regulator problem we shall 

first consider a simple example. Let G, (s) be , and let the 

X (StX) (s+.j) 



1 


I 


I 


r 




Fig. 1.25 


where 

A-Fo l1b*rolc«t2 1]. 

.“3 -4j, [ij. 

This gives us the equivalent formulation of Fig. 1.26, in which the 
initial condition response is considered as an external input. The opera- 
tors in Fig. 1.26 now both satisfy GO = 0. Since it is a first- amd third 

quadrant function, sat kx is a positive operator, and G. (s) is also a 

2 

positive operator, since Re G. (jw) « — + 2 ) ^ 

2 ((D‘^+1) (u +0) 

p(s) + q(s) e s^ + 5s + 5 which is strictly Hurwitz. It is easy to show 
that G^ is strictly positive and bounded, (the finite slope at the origin 
being necessary for this boundedness) . In fact G- is also bounded, 







I 
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Fig. 1.26 

since |Gj^(jw)[ is bounded, since Gj^(s) has no poles on the 

imaginary axis. Now in Fig. 1.26 the initial condition response = 

At 

£ e— Xq is in L 2 because the eigenvalues of A, which are the poles of Gj^(s) , 
lie strictly in the left half plane. We can now apply the Positive Operator 
Theorem as stated above to deduce that the functions Uj^(t) . yj^(t) , U 2 (t) , 
y 2 <t) of Fig. 1.26 all belong to L 2 > This implies that they all approach 
zero as t which is the desired asymptotic stability, provided that 
we can show that all of their derivatives are bounded or square-integrable. 
Now it can be shown for a system of the form 

X « A X + b u 

y « £ X 

that if u e L 2 then y c L 2 , provided that A is an asymptotically stable 
matrix, i.e. all its eigenvalues lie in the half-plane Re s < 0. Thus, 
by describing G^^ in this state space form we see that in Fig. 1.26 y^^ e L 2 . 
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• • 

Then* since C i .< ^2 ^ ^2* operator * sst kx« we 

H ' * * 9 /i ^ ^ 

have ^ G^(x(t)) - ^ • ~ ; thus since <**. 72 ^ L 2 and so Uj^ e L^. 
We therefore ha\^ the desired asymptotic stability. 

Now we consider application of the Positive Operator Theorem 
to the feedback system of Fig. 1.27, which represents our second-order 
lossless voltage regulator. We have replaced the stp x function with a 
sat kx function, k being arbitrarily large, in order to ensure that 
^(sat kx(t)) is bo'inded, as just discussed. 



Now G- (s) «s ■■ ■ " - is not a positive operator, since its Nyquist locus entesrs 
^ s^+1 

the left half plane, as depicted in Fig. 1.20. However, sat kx is 

"very strongly positive ", and by introducing two multiplying factors into 

the loop we can make use of this fact, so that we end up with two positive 

operators. Now G^(s) has poles on the imaginary axis, and therefore its 

initial condition response will not be an L 2 function. To take care of 

this fact we can introduce a very small amount of damping, so that G^(s) 

becomes — for some arbitrarily small r > 0. This represents a small 

s^+rs+1 




1. 
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ii 
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amovint of series resistance in the inductor of Fig. 1.5, which is a case 
we consider in §1.6. 

Let Z{s) be a transfer function, and consider the operator H 
with input u(t) and output y(t) defined by Fig. 1.28. 



Fig. 1.28 


Suppose that 


Then since 


Z(s) 


p(s) 

q(s) 


n . 

qi^s”+...+ qo . 


p(D) u(t) = q(D) w(t) = p(D) y(t) , 
p(D) Iu(t) - y(t) ] = 0 , 

so u(t) = y(t) if euid only if u^^^ (0) ■ y^^^ (0) for 0 ^ i $ n-1. Thus H 

is the identity operator if and only if the initial conditions of its 

second operator Z ^(s) match the initial conditions of the input function. 

B s+l 

Consider now the system of Pig. 1.27 (where G. (s) = ^ '*' ) modi- 

s^+1 

fied to the form of Fig. 1.29, where M is a matrix relating the initial 
conditions of Z ^(s) to those of G^(s ) , as just discussed. We now 
modify this to the form of Fig. 1.30, which still represents the same 
system, provided that Vj^(t) and V 2 (t) aie 0 for all 0 4 t < “, and M is 
suitably chosen. Now let Gj^(s) Z(s) = £(^s~A) so that we have the 
configuration of Fig. 1.31. The configuration of Fig. 1.31 is of the 
form we want, except for the initial condition M x^ in tne feedback opera- 
tor, which we want to be zero so that the feedback operator satisfies GO = 0. 




1 









G^(s) M Z(s) M Z"’(s) 


sat kx 


Fig. 1.29 


G.(s)Z(s) 


sat kx 


Z-' (s) 


Fig. 1.30 


To circumvent this difficulty we assume that the system is started at time 
t = -1 with all initial conditions zero, and that in the interval 0 ^ t $ 1 
the inputs Vj^(t) and V2(t) drive the system to the state shown in Fig. 1.30 
at time t *= 0. We let v^^ and v^ be 0 for 0 ^ t < ". Now the desired 
stability follows from the Positive Operator Theorem as outlined above, 
provided that wc can show that the forward and feedback operators of 
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Fig. 1.31 


Fig. 1.31 are positive, with one of them being strictly positive and 
bounded. We take Z(s) » ys+l where Y ® so that we have 

G (s) Z(s) 

s ^+1 

« 1 + ^ince $Y = 1. 

8*^+1 

i.e. Re Gj^(ju) Z(jw) ® 1 > 0 for all 0); and p(s) + q(s) * 2s^ + (S+Y)s + 2 
which is strictly Hurwitz; thus Gj^(s) Z(s) is a positive real function. 

He now show that the operator of Fig. 1.32 is also positive, where f(o) is 
any first- and third-quadrant nonlinearity with bounded slope at the origin, 
y > 0, and w(0) «» 0. 

Let 

f f(a) da » F(x) , ^ 0 for all x. 

^0 


i 




1 


u =w+yw 


1 

w 

1 +ys 








Fig. 1.32 

F(x) 



Fig. 1.33 


y = f(w) 


fT fT ^ 

I uy dt = I (w+Yw) f(w) dt 
'0 ■'o 

= [ w f(w) dt + y[ f(w) 4r 

rT rw(T) 

= I w f(w) dt + y( 

•'o Jq 

fT 

» w f(w) dt + YF('^(T)) - 
^0 


rT rw(T) 

= I w f(w) dt + y( f(w) dw 

•'o Jq 


f(w) dt + YF(w(T)) - yP(v(0)) 


since we assume w(0) = 0. Thus the operator of Fig. 1.32 is indeed 
positive. 

We make two comments here about the multiplier (1+YS) , Y “ 

Firstly, it is easy to show that in this case this is the only such multi- 

X 

plier (to within a multiplicative constant) which will make — r — become 

8>1 
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positive real. Secondly, we could have predicted that such a multiplier 
exists by making use of a theorem of Brockett and J. L. Willems ([6], 

Theorem 2) , which says that the feedback system of Fig. 1.19 is asymptot- 

* 

ically stable for any feedback gain k satisfying 0 < k < k , if and only 
if there exists a positive real function Z(s) such that Z(s) Ig(s) + ~ 

I k 


is positive real. 

In §1.5 (d) we summarize the application of the Positive Oper- 
ator Theorem to a stability analysis of this kind. 


(e) Stability by Lyapunov's Method 

Now we show that the lossless second-order regulator is stable 
by means of the method of Lyapunov. First, we give some definitions euid 
use them to state the Lyapunov theorem we wish to apply. Then, we make 
use of Dissipative System concepts, as described in Chapter 2, to obtain 
a Lyapunov function for the second-order regulator. Although we have 
already established that the regulator is stable using the Positive Oper- 
ator Theorem, the use of Lyapunov's Method is important here, because of 
its extension to the resistive-source-impedance case, §1.7. A useful 
reference on Lyapunov Theory is [37 ] . 

Consider the system of equations 

x(t) = £(x(t) ,t) ; 0 4 t < 0°; }c(0) = x^j ^(0,t) «» 0^; x(t) e 

A real-valued function V(jc,t) is called a Lyapunov function for this 
system if 

(i) V(jc,t) has first partial derivatives with respect 
to jc and t which are also continuous with respect 
to X and t, and 

(ii) is bounded if | |jc^| | is bc’inded, and 


- . i. 


.j .... 




"T 
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(iii) V(X|t), the time rate of change of V as given by 


"‘■"■(1.%'.) -S 


satisfies 

V(x,t) W(x) ^0 

for all X, and some continuous function WW. 

V(>t,t) is called positive definite if V(0^,t) = 0 and if there exists a 
continuous increasing scalar-valued function of a scalar argument V^(o), 
such that V^(0) - 0 and 

V(x,t) ^ V^(| lx| |) , 

V(5c,t) is called radial ly unbounded if V^ldi approaches infinity as (T 
approaches infinity. V( 2 c#t) is called decrescent if there exists a 
second scalar function of a scalar argument V 2 (CJ) which is continuous 
and nondecreasing r such that V 2 ( 0 ) « 0 and 

V(x,t) ^ V^d |x| |) . 


We then have the following Lyapunov Theorem: 


Theorem 1.3 

If V(5c,t) *= V(x,t+T) is a Lyapunov function for the periodic or time- 
invariant system of equations 

x(t) « « nx,t+T) 

x(t) C 1 r"; 0 5 t < x(0) = = £' 

and if V(>c,t) is positive definite, decrescent, and radially unbounded, 

then the system is globally asymptotically stable about 21 ® provided 

«• 

V(x,t) is not identically zero along any nonzero solution. 


T 


T' 


T 


T 


1 
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For a proof of this theorem the reader is referred to the w^rk of LaSalle 
(17] > and also to [37]. 

We now show that the second-order regulator is stable using 
concepts from the theory of Dissipative Systems, which is closely allied 
with the theory of Positive Operators, as discussed in Chapter 2. 

Let G be an operator (or system) with input u(t) and output 
y(t), 0 < t < “, defined by the equations 



The vector jc(t) is called the state vector for G . Let w(u,y) be a real- 
valued function of u and y. Then G is said to be dissipative with respect 
to the supply-rate w(u,y ) if there exists a nonnegative function VM with 
V(£) = 0 such that 

V - w $ 0 . 

V is called the storage function for G . In this thesis the only supply 
rate w which we consider is the product uy, so that when we refer to an 
operator as dissipative we have this supply rate in mind. In Chapter 2 
we shall show that an operator is dissipative if and only if it is a 
positive operator. 

Consider for instemce a one-port KLC electrical network with 
input current u(t) and voltage y(t), as in Fig. 1.34. The operator G 
mapping u(t) into y(t) represents the impedance of the network; expressed 
as a transfer function G(s) it is a positive real function. The requirement 



y(t)f 


RLC 

network 


Fig. 1.34 


for G to be a positive operator 


•-C 


uy dt 5 says that the net 


energy absorbed by the network is ut all times positive, the concurrent 
requirement being that G maps zero into zero, i.e. that nu initial energy 
is stored in the network. If G is viewed as a dissipative operator the 
storage function V is the energy stored in the network, and the supply 
rate w •« uy is the rate at which energy is supplied. The requirement 
V - uy 5 0 says that the rate of increase of stored energy is not greater 
than the rate of supply, (because of dissipation within the network) . This 
requirement is independent of the initial conditions within the network. 

Now consider the feedback connectibn of Fig. 1.21 with 
Vj^ = Vj = 0. Let Gj^be a dissipative operator with state vector and 
storage function G., be dissipative with amd ^ 2 ^— 2 ^ ' 


Vj - s 0 


Vj - ujyj « 0 


Uj = au.d '*2 “ ^2' " '^1^—1^ ^ ^2^—2^ 


obtain 


V - . Vj 


- Vj - . Vj - Ujyj 






\ 


Thus* a Positive Operator Theorem proof of stability provides us with a 


Lyapunov function for use with the Lyapunov Theorem, provided that the 
positive operators can be described in dissipative operator form. 

Fig. 1.35 rep'‘csents the second-order regulator, with 0 • 1 
for notational convenience. 



Fig. 1.35 


A state-space description of this system is 


0 -1 


1 0 


+ 1 


II 11 I X 


-stp z 


where 


X <■ A X b u 


z « £ JC 


u ■ -f(z) 


£(IS-A) 




1 


? 


1 


1 


Fig. 1.36 depicts the same system with multiplier (l-i-s) included in the 
forward path and (l+s) ^ Included in the feedback path. As we saw in 
§1.4 (d ) , introducing these factors has no effect provided that the 
initial conditions of /-rr^ are correctly chosen. This means that ther« 


fe) 


and those of (1 + 


cf correctly chosen. This means that there 

is a linear relationship between the initial conditions of the operator 
/*' 



Fig. 1.36 

A 8tate*space description of the system of Fig. 1.36 is a non*minimal 
realization of the system of Fig. 1.35. We have 


0-1 + 1 u 

1 0 0 


[2 01 I *1 I + 


-z + y with z(0) • Xj^(O) + XgCO) 
-stp z 


i.e. 



where 


Note that 


since 


X “ A X + b u 


y • £ X d u 


z « -X + y 


u ■ -f(z) 


cds-a)”^ b + d ■ 1 + 

“ *■ s^+1 


z(t) = Xj(t) + Xj(t) , 


'z + z ■ y ■ 2xj + u » (Xj^+X 2 > + (Xj^+Xj) . 


From Chapter 2 we know that a storage function V. (x) for (1 + — ) is 
1 

— £'K }c where K ** K* satisfies the matrix equation 

b.‘b * (2d)“^(b'iC - c) ■ £ . 

In this case this equation has a unique solution iC •* 2^, so that the storage 
1 + j is Vj^(x) «= x'x ■ Xj^ + Xj • In Chapter 2 we show 
that the urdque storage fxinction for the feedback operator of Fig. 1.37 is 
sod z ■ z stp z. 



Fig. 1*37 


Since z = Xj^ + Xj, we have VjCsc) = Sod(Xj^+X 2 ) • Thus a Lyapunov function 
for the second-order regulator described by 
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Xj^ « -X 2 - stp(Xj^+X2) 
^2 “ 


is 


V(x) = V^(x) + V^(x) « x^^ + x^^ + Sod(x^+X 2 > 
Differentiating we get 

VW = + SXjX^ + (Xj^+x^) stpCXj^+x^) 

« - 2Xj^ stp(x^+X2) + 2 Xj^X 2 + (Xj^-X2-stp(Xj^+X2) ) stp(Xj^+x,) 

2 

«= “(Xj^+X^) Stp(Xj^4x2) - Stp (Xj^+Xj) 

2 

» - SodCXj^+x^) - stp (Xj^+Xj) 

<0, *» 0 if and only if x^+x^*©. 

But 

• • 

Xj^ + x^ = - X 2 - stp(Xj^+X2) 

« 0 on {Xj^+X 2 = 0 } only at Xj = X 2 = 0 , 

i.e. V(x) is identically zero only along the solution x “ Thus, by 
the Lyapunov Theorem 1.3 we obtain the desired global asymptotic stability 
about 2L “ — * Notice that for Fig. 1.36 we have proved b-ability for any 
initial conditions in the forward path and any initial conditions in the 
feedback path; this includes the case where there is a linear relationship 
between these two sets of initial conditions, this case being the system 
of Fxg. 1.35. 


In conclusion of this section we note that the Lyapunov method 
is superior to the Positive Operator Theorem method for determining 
stability, since with the Lyapunov method initial conditions do not require 
special attention, nor do we have to modify Gj^(s) to have poles which are 
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Off the imaginary axis. The Positive Operator Theorem is perhaps more 
useful for providing a Lyapunov function. We also found it necessary to 
modify stp x to be continuous at the origin, however the behavior at the 
origin is taken care of by the chattering behavior analysis of §1.4 (b) , 
in which stp is better left as a discontinuous function. The Positive 
Operator Theorem and the Lyapunov method are thus used to deterrrlr.e that 
the state will reach a neighborhood of the origin from any starting point 
in state space. 

§1.5 Fourth-Order Lossless Regulator 
(a) Choice of a Feedback Law 

The second member of the series of voltage regulators of Figs. 
1.1 and 1.2 is shown in Fig. 1.38, where we have assumed unit values for 
the components. As with the second-order regulator this assvur.ption barely 
limits the generality of our case; it represents all regulators of the 
type of Figs. 1.1 and 1.2 where the inductances are all equal and the 
capacitances are all equal. We want lim z. (t) « a for some given 

t -► 00 

0 < a < 1. 
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The state evolution equations are 


r • 

*1 


‘ 0 

-1 

0 

0 ' 


'*1' 

+ 

' 1 ' 

• 

*2 


1 

0 

-1 
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*2 


0 

• 

*3 


0 

1 

0 

-1 


*3 


0 

• N 
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0 

0 

1 
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_ 0 


• • 

Now in steady-state operation we want = a, and Z2 = = 0. The 

natural ext''nsion of the control law we chose in §1.4 (a) is thus 

u = (1 , if $jZ^ + ^ ° 

( 0, otherwise, 

where 62 » are suitably chosen positive constants. 

(b) S tability by Total Gain Linearization 

We make the change of variables = 2z^ , X2 = Hz^o.) , 

“ 2Zj , x^ = 2(z^-a) , and thus obtain 


*1 

= 

“ 0 

-1 

0 

0 ' 



+ 

" 2 { u - a ) ■ 

*2 


1 

0 

-1 

0 


*2 


0 

^3 


0 

1 

0 

-1 


^3 


0 

*4 


0 

0 

1 

0 


*4 


0 


where 2(u-a) = -stp (B,x + B x + B^x, + B .x.) . By expressing 

(-2+2a,2a) ^ 

Xj^, x^, x^ in terms of x^ = x we can put this in the form 

‘x*+ 3x + X = -stptBj^ X* + B2X + ^B 2 +B^)x] . 

The desired stability is about x - 0. Taking Laplace transforms we obtain 
the feedback representation of Fig. 1.39. This is of the form of Fig. 1.16. 
Using the Routh-Hurwitz Criterion it is easy to show that the associated 
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feedback system of Fig. 1.19 is stable for all 0 < k < “ if and only if 


S3 S 2S3 


63 <63(63*63) 


6364 « 6363 . 

In particular we note that 3=3^ = 6.. = 34, = 1 satisfies these require- 

1234 

ments. In §1.5 (d) we shall show that indeed Aizerman's Conjv,c:ture is 

correct here for 3^^ = 82 - 3^ = 3^ “ 1. This may not be too surprising, 

s^+s^+2s+2 

considering the nature of the Nyquist Locus of G(s) = — , which is 

s>3s^+l 

shown in Fig. 1.40. The pole-zero pattern of G(s) is shown in Fig. 1,41. 
For design purposes the total gain linearization method appears to provide 
a reasonable approach for an (initial) investigation of the stability of 
this class of regulators; the above three inequalities tor 3w 6 ^, 3-,# 3 „ 
arc probably the necessary and sufficient conditions for global asymptotic 
stability of the fourth-order regulator. 


I 




...A .. 






Fig. 1.40 


jcu 



Pig. 1.41 

(c) Chattering Behavior 

We follow the development of § 1.4(b). The system of Fig. 1.39 
is of the form of that in Pig. 1.15, with 
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-1 0 O b 


0 -X 0 


1 c 


[®1 ^2 ^ 3“®1 ^ 4 ] 


10-1 


0 1 0 , 


a «s 2 - 20 , b = 2a . 


The switching surface S is the hyperplane £ x “ » 


+ $ 2*2 = 0 , 


and X C S is an endpoint if 


-acb<cAx<bcb 


-(2-20)Bj^< $2*1 + (33-2 Bj^)x2 + *3 * <2l“®3)*4 < • 

For example when = 62 “ 3^ = 6^ = 1 we have: S is the surface 

Xi + X2 + X4 = 0, and the conditions for chattering to occur at ac £ S are 

-2 + 20 < Xj^ - X 2 < 20 


-2 + 2o < -2x. - X. < 20 
2 4 


2o < 2Z2 + Z4 < 20 + 1 


Motion along S is governed by ac = F x where 


detUs-P) ■ |“t33S^ + 628^ + (B3+B3 )s + (32+34)! 


3 2 

Now the conditions for a cubic polynomial a^s + a2S + aj^s + a^ to be 
strictly Hurwitz are 




} 


! 


f 


*0’ *1' *2' “3’ ° 


* 2*1 - * 3"0 ’ “ • 


Therefore for asymptotic stability of the chattering mode on S we require 


61, 63, 6^ > 0 


S263 > B36, 


The second of these is a strict inequality (whereas in (b) above we had 

and is not quite satisfied by 6, = = 6-, = *= is for 

1234 

these values of we have 

det(Is-F) - s(s^+s^+2s+2) 


s (s+1) (s +1) . 


However in what follows we use these borderline values for convenience. 

(d) Stability by the Positive Operator Theorem 

Following the development of §1.4 (d) , we wish to find a multi- 
plier Z(s) such that G(s)Z(s) is positive real, where 


G(s) = 


3 2 

s +s +2S+2 


(s-H) (s +2) 


s^+3s^+l (s^+2.61) (s^+0.i8) 


We also want Z (s) followed by stp x to yield a positive operator. At 
this stage we encounter the difficulty of testing a function for positive 
realness. There are several different characterizations of positive real 
functions, ([16] Chapter 5), but it seems that the least difficult method 
to apply is to test whether Re G(j(jj) > 0 for all u, and p(s) + q(s) is 


strictly Hurwitz. Using this method one can show thft for 
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1 


1 


1 


8^+s^+2s+2 

G(s) ■ - ^ the required multiplier 2(s) cannot be of the form (s+a ) , 

8 +3s^+l 

nor of the form - ■. The simplest multiplier Z(s) is therefore of 

3 2 

, s +as +bs+c _ ... ... 

the form ' ' ■« For this we obtain 

s^'+ds+e 


Re G(jw) Z(jw) = 


(2-0) ) H(fa)) 

(W^-2.61) (w^-0.38) + d^to^] 


where 


H(u» « (c-aw^) (e+dw^-u)^) + (b-ts^ (d-e+U)^) . 


Thus Re G(joj) Z(jw) J 0 for all w if and only if 

H(w) « k(2-w^) (w^-2.61) (w^-0.38) for some k > 0 , 


and this leads to the requirements 


c e = 2 


b(d-e) + c(d-l) - a e = -7 


b+e-d-ad+a^S 


To ensure that Z (s) followed by stp x yields a positive operator we 
make use of the following theorem of O'Shea {[26], [32], [12]); 


Theorem 1.4 

The operator F shown in Fig. 1.42 with input u(t) and output y(t) and Z(s) 
rational, is a positive operator for any monotone nonlinearity f(a) for 
which f(0) =0, if and only if 


z(s) “ g^ + Ys - g(s) 


where 


Y ^ 0 , g(s) = g(t)e dt , g(t) ^ 0 , and g ^ g(0) * g(t) dt. 
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Fig. 1.42 


For a proof of this theorem the reader is referred to reference [12] • 
Applying this theorem to our Z(s) we obtain the constraints: 

| a> b, c, d, e > 0 “(4) 

d ^ a “(5) 

d^ ^ 4e -(6) 

d(a-d) > b-e -(7) 

e(a-d) > Y "(8) • 

After some trial and error one can find a set of values for a, b, c, d, e 
which satisfy (1) through (8) . One such set of values is 

a ■ 200 ,b=96,c®4,d* 1.45 , e = 0.5; 

so that 

. s^+200s^+96s+4 

z(s) 5 . 

s +1.45S+0.5 

We also need to show that G(s) 2(s) » has p(s)+q(s) strictly Hurwitz. 

p(s) 

By following through an argument similar to that of §1.4 (d) we conclude 
that the fourth-order regulator of §1.5 (a) is globally asymptotically 
stable when B ® ® 6.. ® 6- = 1. The practical implications of this are 
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that it is just as easy to stabilize a fourth-order regulator as it is to 
stabilize a second-order regulator. 

To sununarize the application of the Positive Operator Theorem 
to switched regulators of this kind, we see that the first step is to 
bring the problem to the fonit of determining stability of the null state 
for a system of the type shown in Fig. 1.15 or Fig. 1.18. We then consider 
the associated system of Fig. 1.19. If this is not stable for all k > 0 
then we cannot make any conclusions in general, though instability for all 
k > ko > 0 will imply local instability in Fig. 1.15, and instability for 
all k < will imply instability in Fig. 1.15 for large initial conditions. 

If the system of Fig. *1.19 is stable for all k > 0 then we know (10 ], 

Tlicorem 2) that tliere exists a class of positive real functions ^ such that 
G(s) Z(s) is positive real for each Z(s) in Stability then follows if 
there is a Z^(s) in ^ such that 2j^(s) satisfies the conditions of the 0*Shca 
Theorem alcove. Behavior in the chattering mode is analyzed as outlined in 
§1.4 (b). 

(e) Stability by Lyapunovas Method 

In §1.4 (e) we saw how to obtain a Lyapunov function for the 
feedback system, by making use of the fact that we could express the forward 
and feedback positive operators in dissipative form. For the linear operator 
given by G(s) Z(s) = £(XS“^) ^ b + d we need to solve the matrix equation 

Ji ii + A'ii ^ k ■ £') (2d)"^(b'K-c) •= 0 . 

This can usually be done, using suitable numerical methods if necessary. 
However, to obtain a storage function V(x) for an operator of the form of 
Fig. 1.42 is not easy. In Chapter 2 wc shall show that if a positive 
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operator maps u(t) into y(t) , the functions 


V (X) 

a — 


sup ~ f W(t) dt , V (x) » inf /I uy dt 

•'o ' ■ 

x(0)“x ^ 


- “ x(0)*0 

are suitable storage functions. Unfortunately these definitions 
require the solution of a nonlinear optimization problem which can be 
solved for all initial conditions only when Z(s) is of order 1 or 2. For 
our fourth-order regulator problem Z(s) is third order, so we mush try 
some other means to find a suitable V(}c) . The most promising approach 
!5eems to be to attempt to find a realization (A, b, c) of Z(s) which 
satisfies the requirements of the following theorem: 


Tl ieorem 1.5 

Suppose -A is a hyperdominant matrix, b' ■ [0 0 ... 0 X] , X > 0 , 

£ = to 0 ... 0 1], and f is any monotone nonlinear xty with f(0) « 0. 

Then the operator mapping u(t) into y(t) defined by 


X “ A X + b u 


z “ c X 


y = f(z) 


1 r 

is dissipative, with storage function V(x) “ T 4 F(x.), where 

i-1 ^ 


f(o) da . 


We shall give a proof cf this Theorem in Chapter 2. A hyperdominant 
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Fig. 1.43 

n 

matrix M is one for which m. . $ 0 whenever i ^ •j, and ^ m. . ^ 0 and 
n ” i-1 ^ 

I m. . 5 0 for all l»j‘ 
j*l 

For our fourth-order regulator we have 


z“^(s) « H(s) 


s -I- 1.45s 0.5 

+ 200s^ + 96s + 4 


The standard controllc±ile realization ([" ] Chapter 17) is H(s) ■ htis-P) £ 
where 

F=ro 1 ol£*ro1 h to.5 1.45 1] 

0 0 1 0 
-4 -96 -200 , 1 , 


Now let 


1 0 
0 1 


p"^ - 


[0.5 1.45 IJ [-0.5 -1.45 1 

so tljat we obtain the realization H(s)= £(l.s-A) ^ where 
. ^ — 1 r*.. . — PaT — 


0 1 0 
-0.5 -1.45 1 

95.3 192.4 -198.6 , 


b “ P^ £ 


c-hP*-I0 0 1] 


L 
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Now -A is not hyperdominant, but if we choose R such that R b ■ b and 
•1 “1 

£ R ■ £ then ~R A R may be hyperdominant. Such an R is of the form 


b o') 


d -b 


ad “ be ■ 1 


c d 0 


0 0 11 . 


Xt seems to be an inpossible task to find values for a. b. c. d which will 
make “R A R ^ hyperdominant, without the aid of a digital computer. In 
Chapter 2 we reconsider this question using tridiagonal realizations- We 

also show in Theorem 2.7 that it would even be helpful to find a, b, c, d 

-1 ” 
which will make -F. A R M column dominant, l.e. m.. ^ T |in. .1. 

- jj ' iji 

i»<j 

(f) Higher-Order Regulators 

As a design procedure for sixth- and hi'<sr-order regulators 
it seems reasonable to assume that Aizerman's Conjecture holds true for 
the class of feedback systems obtained by using the kind of feedback control 
law described in §1.4 (a) and §1.5 (a) . The chattering mode analysis of 
§1.4 (b) and §1.5 (c) shows that the ntimerator polynomial of the forward- 
path transfer function must be strict Hurwitz, but this reqtiirement will 
be covered by the conditions obtained by using total gain linearization, 
(i.e. the system of Fig. 1.19 must be stable for all 0 < k < *) . 


§1.6 Second-Order Regulator with Inductor Loss 

In this section we consider the regulator of Fig. 1.5 with a resistance 
added in series with the inductor. We know in advance that this will not 
affect our conclusions about the stability of the second-order regulator, 
since "Dissipation aids stabilization". However, we wish to obtain a 


Lyapunov function for this case, in preparation for §1.7. 
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Fig. 1.44 shews the regulator under consideraticn« with unit values 
for source voltage » inductance* and capacitance. The series resistance 
has value r. 



Fig. 1.44 


T^ollowing the pattern of §1.4 and §1.5 we obtain the feedback representation 
of Fig. 1.45* which is seen to be of the form of Fig. 1.18. As in §1.4 (e) 



Fig. 1.45 

we take B«»l for notational convenience. Foi r > 1, G(s) ■ is a 

s^+rs+1 

positive real function* so that no multiplier is necessary in such a 
case to prove stability using the Positive Operator Theorem. Fig. 1.46 
depicts the Nyquist locus of G(s) when ^ ® ' 2 » However we do not 

make use of this positivity of G(s) for r ^ 2* because we want one storage 
function VU) for all r* for the purposes of §1.7. We must use the same 




Fig. 1.46 


(s+ir 


multiplier as before, Z(s) = s + 1. Then G(s) Z(s) = 

s +rs+l 


, which is 


a positive real function for all r ^ 0. We have 


G(s) Z(s) = 1 + = c(Is-A)"S 3 + d 


s +rs+l 


where 


A - 


-r -1 
1 0 


1 

0 


c “ [2-r 0] , d = 1 


As in §1.4 (e) we therefore describe the system of Fig. 1.45 by 

X e’ A X + b u 


y “ £ X + du 


z = -z + y with z(0) = Xj^(O) + X 2 ( 0 ) 
u = -stp z 


In Chapter 2 we introduce the concept of a storage function and show that 
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a storage function V, (x) for (1 + is given by t x'K x where K ■ K' 

” s +1 ^ 

satisfies the matrix equation 

K A + A'K + (K b - o') (2d)‘^(b'K-c) = 0 . 


If we let K 


ki 

*^2 *^3 


tiien we obtain 


k-‘ - 4k, = 0 

A 4» 


k2(ki+r-2) + 2(k3~rk2-kj^) « 0 


(kj^+r-2)‘ + 4(k2-rk^) « 0 


from which we get the four possible solutions 


= r r+2-2/2r 


r+2-2v^ 


= [~ r+2+2>^ 


r+2+2/2r 


Kj * r+2-2/2r-4 


r+2+2/2r-4 


^ = r r+2+2/2r-4 


r+2-2/2r-4 


We see that for r < 2 the solutions K« and K, are not real# so we consider 

—3 

only and K^. Since we see that j must be the available 

storage function ^ and ^ 21* ^^ 2 — ^^^st be the required supply (see Chapter 2) . 
Since the convex combination of these two storage functions is again a 
storage function, we know that ^ x’2. 21 ^ storage function whenever 

2 . “ (l-n)K 2 for some 0 $ n $ !• Thus, if -1 ^ Y $ 1 ve have 

the storage function 

V^(x) » i (r+2+2y/2‘r) (xj+x^) . 


To check that the operator mapping u into y is dissipative we evaluate 
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V - uy ■ -j(r+2+2Y»^) (Xj^Xj^+XjX^) - (2-r)Xj^u - u^ 

■ - t(r+Yv^)Xj^ - u]^ - (l-Y^)Xj^^ on substituting 
$ 0 since |y| ^ 1. 

We can now establish stability by noting that the storage function for the 
feedback operator is as it was in §1.4 (e) , namely V^Cjc) * Sod(Xj^+X 2 ) • 

Thus if we describe the regulator of Figs. 1.44 and 1.45 by the equations 

■ 

= - r Xj^ - x^ - stpCxj^+x^) 

*2 “ *1 

we have the Lyapunov function 

V(x) « Vj^(x) + V^U) = i(r+2+2Y>^) (Xj^'+Xj^) + SodCXj^+x^) . 

Then 

V(x) = (r+2+2Yv'2r) (x^Xj^+x^x^) + (x^^+x^) stp(Xj^+X 2 ) 

e - t(r+Y»^)Xj^ + stpCXj^+x^) - 2r(l-Y^)Xj^^ - SodCx^^+x^) . 

By choosing Y -1 we have 

V(jt^) $ 0 , = 0 only at Xj^ " x^ 0 . 

This allows us to conclude stability about x = £ \ising the following theorem 
of Yoshizawa P9 j , whiuh we introduce hers with a view to the time-varying 
situation of §1.8. 

Theorem 1.6 

If V(x,t) is a Lyapunov function (as defined in §1.4 (e)) for the system 
of equations 
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x(t) “ 

x(t) e IR*' ; 0 ^ t < “ j x(0) « x^ } f(0,t) * 0 , 

and if V(x#t) is positive definite, decrescent, and radially xinbounded, 
then the system is globally asymptotically stable about x * provided 
that -W(x) is positive definite, where V(x,t) < W(x) < 0. 

For a proof of Theorem 1.6 the reader is referred to [37]. 

1.7 Second-Order Regulator with Resistive Source Impedance 
(a) Introductory 

Inclusion of a resistance in series with the source voltage E 
of Fig. 1.1 leads to a much more difficult stability analysis. For the 
second-order regulator with source resistance, the network acts some of 
the time (u = 0) like that of Pig. 1.5, and for the rest of the time 
(u « 1) like that of Fig. 1.44. Hence, using the same control law we 
might expect to have stability, since both of the regulators of Figs. 1.5 
and 1.44 are stable. Furthermore, stability is to be expected from the 
notion "Dissipation aids stabilization". In fact, we do find that these 
expectations are true; however, this not so easy to establish. For 
the second-order case we can fall back jn a phase plane analysis, but for 
higher-order cases a Lyapunov approach seems the only way. We provide in 
part (c) here a Lyapunov analysis of the second-order case, which makes 
an elegant use of the ideas of dissipative systems. In part (b) we apply 
a total-gain linearization as an initial investigation of stability. 

The second-order regulator is shown in Fig. 1.47. 








If II ■ fulfill^ nil , ■■iiiAoaiii^ 








T 





I 


l^u 



1-*-% 

T -*2 


rig. 1.47 


The evolution equations are 


0 -1 I + u| -r 


■' ° 1 I [ 

0 0 J I [ z. 


♦ I u 


We see that the control variable u enters the state evolution equations 
in a multiplicative way, as well as in the usual additive way. As before 


1 , Zj + 2^ < a 
0 , 22 + 22 > a 


and X, 


1 ' 2 


2 z 2 “ 2a, which yields 


x^ “ -X2 - urXj^ + 2(u-a) 


X2 = 


2(u-a) ■ -stp 

(-2+20, 2a) 




j 
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» "2aqx^ " ^2 ^ (qx^-1) stpCXj^+x^) 
*2 “ 


where ^ notational convenience. Now we can write this equation in 

two ways^ the first way being 

*. • • • 

X + 2aqx + X = (qx“l) stp(x+x) 


and the second way being 

•* * * * 

x + ql2a “ £tp(x+x)l X + X “ - stp(x+x) 


•• • • 

X + rfjiW X X = - stp(x+x) 

where <t)(}c) = -^t2a-stp (x.+x )]. i})(x) takes on the values 0 and 1. 

^ (-2+20, 2a) ^ ^ - 

The first way leads us to the feedback system representation 
of Fig. 1.48. On comparing this with that of Fig. 1.45 we might hope that 
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the feedback operator in Fig. 1.48 is a positive operator. However, this 
is not the case. The representation of Fig. 1.48 does not seem to be 
helpful at all, so we look at the other way of describing this regulator. 
This second way in state-space form is 



Fig. 1.49 depicts this in feedback form: 



Fig. 1.49 





-70- 


(b) Stability by Total Gain Linearization 

Following the total-gain linearization argument of §1.4 (cj , 
and comparing F^g. 1.49 with Fig. 1.17, we expect stability for the 
system of Fig. 1.49 if it is stable when stp 0 is replaced by ko, for all 
0 < k < “. This gives us the equations 


*1 

= 

~-(2aq+k) 

-d+k) " 


"^l‘ 

+ 

qkXj^(x^+X2> 

*2- 


1 

0 


-*2- 


0 


X 


A 3C + ^(jc) 


For these we can use the following theorem: 


Theorem 1.7 

The system of equations 

3i(t) = A(t) x(t) + f,(x,t) 

x(t) C m" ; 0 ^ t < « ; x(0) = 3^ ; £(0,t) = 0 

is locally asymptotically stable (i.e. if IIjCqM ^ M for some 0 < M < ®) if 

X(t) = A(t) X(t) , X(0) = JCq 

is exponentially stable, and 

Lim ■■■ II ,1 = 0 , uniformly in t. 

||x|| - 0 ilili 

By Mill mean the usual Euclidean norm exponential stability means 

that 1 lx(t) II < ae for some aib > 0; and by 

Lim g(x,y) = 0 uniformly in y 
0 


1 


I 


I 


! 


we mean that for all h > 0 there Is a £ > 0 which is independenc of y, such 

that |g(x,y) | < n whenever jx] <6. For a proof of Theorem 1.7 the reader 

is referred to Chapter 4 of Bellman [ 4 1 . 

2 

Now we have det(^s-A) ■ s + (2aq+k)s + (1+k) , thus A has its 
eigenvalues in the left half plane, so 21 “ A )C is exponentially stable. And 


II£(2.) I 1 jqkx^(x^+X2)| 


0 as X 0, for all 0 < k < 


/ 2 . 2 

^Xi +X2 


Thus we have local asymptotic stability for the system when stp a is 
replaced by ko. 

(c) Stability by Lyapunov's Method 

In attempting to find a suitable Lyapunov function we can start 
with that of §1.4 (e) and modify it by trial and error. Here let us take 
a a: for convenience, so that stp becomes sgn. In §1.4 (e) we had 

V()c) » — (x +x ) + — |x+x|. One possible approach here is to add to this 
a term in |x|. So let us try 

V(x) ■ ■j(x^+x^) + a|x+xl 4 b|x| where a,b > 0. 


• • ••• 


V = XX 4 XX 4 a(x4x) sgn(x4x) 4 bx sgn x 


• 2 • • • 2 • 

-qx [l-sgn(x4x)] 4 [(a-l)x-ax] sgn(x4x) - a sgn (x4x) 


■ 2 • 

4 aqxtsgn (x4x) - sgn(x4x)] 4 bx sgn x . 


Now if we take a ® “ 2aq we have 

• *2 • li*il 2**1* 2* • 

V «s -qx [l-sgn(x4x)l - — |x4x| - sgn (x4x) 4 ^xlsgn (x4x)-sgn(x4x) 42sgnx] 

1 I * I 1 * * * 

^ W «= - — |x4x| 4 •- qx (l"Sgn(x4x) 4 2 sgn x] when x 4 x 0. 


1 


! 


1 


1 
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> • > 

By looking at the four cases given by x < x ^ 0, we find W ^ 0 (and 

• * 1 

hence V ^ 0) for all x#x if and only if q < ~. Thus we have the desired 

global asymptotic stability, using Theorem 1.6, provided that r < 1. This 
makes us wonder whether we might obtain instability if r is large enough. 

Further study of the above V shows that we can improve on 
this bound on r. For, collecting terms in a different way wo obtain 

• «2 • • * 2 * 

V = -qx tl-sgn(x+x)] + ( (a-aq-l)x-ax] sgn(x+x) - a sgn (x+x) 

® 2 

+ aqx sgn (x+x) + bx sgn x . 

Now if we let b » aq and a-aq-1 * -a, i.e., a « 2^' requires q < 2 

for a > 0) , we get 

V = -qx Il-sgn(x+x] - a|x+x| - a sgn (x+x) + aqx[sgn (x+x) + sgn x] 

1*1 • 2 • 

^ W B - a|x+x| + aqx (sgn (x+x) + sgn x] . 

> • > 

Looking at the four cases x < 0, x < 0, we now find that W ^ 0 (and hence 
V < 0) for all x,x i^ only if q < 1, so by Theorem 1.6 we now have 
global asymptotic stability provided r < 2. 

We can do better than this, however, using the ideas of 
Dissipative Systems. We make use of the fact that for a dissipative 
system, the convex combination of V and V is also a storage function, 
where V is the available storage and V is the required supply (see 
Chapter 2)* Now in §1.6 we described the introduction of the multiplier 
(s+1) by a nonxninimal state-space representation. If we follow that 
method hero we obtain the nonininimal state-space representation: 

I 




«,i 


T 


T 


T 


T 
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1 


1 


1 


is a dissipative operator. Let us try the storage function we obtained 


in §X.6^ i.e. 


V(x) » ■^(r+2+2y*^) (Xj^ +Xj ) where $ Y $ 1 


V “ uy «» (r+2+2Y*^) (x^Xj^+X 2 X 2 > ” uy 

« -r(|»(r+2+2Yv'2r) + [r(l+4)) + 2Y*^1 x^^u - u^ 

I -t(r+Y*^) Xj^ - u]^ - (1-Y^) Xj^^ when 4> - 1 
(r+2Y»^) Xj^u - u when 4i - 0 


We will therefore always have V - uy ^ 0 if we can choose Y ®° that 
(r+2Y*^) “ 0, i.e., Y ■ “ Since Y ^ “1 this means that our ^.pcraLor 
is dissipative for all r 5 8. The storage function for the forward-path 


2 . 2 


operator is therefore = (x^^ +x^ ) , and for the feedback operator of 

Fig. 1.50 we know the storage function is •* Sod z “ Sod(Xj^+X 2 ) . Thus 

2 2 

the Lyapunov function to try is V “ (Xj^ +X 2 ) t Sod(Xj^+X2> , which, admittedly, 
is the one we started with. Indeed, since 


we have 


X + r(Mx) X t X ■ - stp(x+x) 


V*=2x X+2XX + (x+x) stp(x+x) 

«s -2r((ix^ - [(r(j)+l) x + x] stp(x+x) - stp^(x+x) 

■ -2r(}ix^ - r(})X stp(x+x) - stp^(x+x) - Sod (x+x) 

■ - X + stp(x+x)]^ - ^(8-r4>) x^ - Sod(x+x) 
$ 0 for all x,x if r $ 8, since 4> i® 0 or 1. 
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When ♦ • 0 , V ■ 0 x + x * 0 , so we need to use Theorem 1 . 3 to conclude 
stability* which we have since x x ■ 0 is not a trajectory of the system. 

Me have not been able to obtain a Lyapunov function which gives 
stability for all r ^ 0. Perhaps the next step in this direction is to 
try to find such a V by looking directly at the definition of a storage 
function. 


(d) Stability by Phase Plane Analysis 

Fig. 1.51 shews a set of phass-plane trajectories for the case 
•j, r ■ 10. l.e.f the trajectories of Fig. 1.51 are the solutions of 


X + 10[---- ^ 1 X + X ■ -sgn(x+x) . 


In order to obtain Pig. 1.51 we need to know the shape of the trajectories 
for 

X + rx + X = 0 . 


These are shown in Fig. 1.52, for r « 10. For a discussion of phase-plane 
trajectories of this type, the reader is referred to Chapter 7 of refer- 
ence [11] and Chapter 3 of reference [37] . 

Prom Fig. 1.51 we conclude that all trajectories will reach 
the chattering region of the switching line. To determine whether the 
resulting chattering motion is stable we need to apply the methods 
developed in §1.4 (b) . Our system obeys 

U ■ - stp £ £ 

(-a,b) 




where 


1 


1 
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A(x) 


-r<Mic) -1 

1 0 


♦ (x) - <|-[2o - stp <Bx,+x-)] 

(-2+2a,2a) 


1 

0 


£ “ 18 1] a » 2 - 2a b *s ?a 


Let 




s 

r 0 -1 1 

and A, <= 

‘ -r -1 ' 



—1 



1 — 

M 

O 

1 


1 0 


Then 


£ “ ^ 3C - b b and £ “ iL ® ^ 


The switching line S is £ x = 0, i.e. 8xj^ + x^ * 0. The first endpoint 
condition is £ <0, which on substitution becomes 

X 

2 «2. 


8%1 


The second endpoint condition is £ £ <0, which becomes 

< e-?°l i! i£ r < 2±i 

2 8+l-r8 8 ' 

> -S 2-2a) jJ, if r > 6±i 

Tl.uSf for Cl = 8 = 1» r > 2, as in Fig. 1.51, we have diattering behavior 


on the half line 


> max(- 


Motion along the chattering line obeys 


X X c f 

X = f - (f"-f ) 


cCf'-f"^) 


i 


1 


J 



I 


1 


I 


f 


which « on substituting, becomes 


“ -(^) , with + 

-*oJ L*,. 


X 2 “ 0. 


Thir is the same result as for §1.4 (b) . We therefore have the desired 
global as^nraptotic stability, for all B > 0 and for all r ^ 0. 

§1*8 Further Refinements 
(a) Preregulation 

The source imp. 'dance shown in Fig. 1.47 is the simplest kind. 
It does not increase the dime*ision of the state equations. In practice 
the source impedance may well be like that depicted in Fig. 1.53. The 
capacitor C provides a form of preregulation , that is, a smoothing of ary 

1 _ 11 L I4 





- v,:. 

nC 





Fig. 1.53 

time-Veu:iations in the source voltage. In that case the evolution equations 




= 2''2 - 


CjVj = G(E-Vj) - ul 


1 


1 


1 



Suppose, for example, that to^ ■ “ b « 1, and that we weuit stability 
about Zj ® Ct. Then, letting = 2Zj^, x^ * 2 X 2 - 2a, x^ * 2z^ - 2, we 
obtain 



Let 2u - 2a * -stp (Sx.+x.) , as before. Then we have 

(-2+2a,2a) ^ 



where 

•{•(x) = •|•t2a - stp (3x +x )] 

(-2+2a,2a) 

We now wish to prove that this defines a motion which is asymptotically 
stable about x, = 0. As yet, we have net done this. Probably the desired 
stability can be obtained by extending the methods of §1.7 (b) and (c) , 
that is, to view these equations in the form of a system like that of 
Fig. 1.49. Certainly it is to be expected tha* this system will be 
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Stable. An even more difficult stability question will be provided by 
replacing the resistor of Fig. 1.53 with an inductor. 

(b) Time-Varying Source Voltage 

So far we have only considered a fixed source voltage, and have 
shewn that the output voltage will settle to its desired value after a 
change in input voltage or a change in the load resistance. In practice 
the source voltage will be time- varying, between limits. Usually this 
variation will be periodic, as for instemce in rectification applications. 
In terms of the regulator of Fig. 1.5, this means that we have 

E(t) ■ Eq + Ej^(t) , where E^ > 0 

and Ej^(t) ^ 0 for all t. 

This modifies the evolution equations of §1.4 (a) for and to be 


~K~ 

cr 

■ 0 -1 ■ 


' *l“ 

+ 

‘ l+e(t)‘ 

-*2- 


_ 1 0 


-^2 


0 


where e(t) ^ 0 for all t. Using the same feedback law we obtain, for 


and x^ 


i.e. 


where 


X = -X- - stp (8x +X-) + e(t) [2a - stp (Bx.+x )] 

^ ^ (-2+2a,2a) ^ (-2+2a,2a) 


^2 “ *1 



f (£ X,t) 

0 


i 


] 


i 


i 
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f(c x,t) 


-2(l+e(t)) + 2a , •« < 0 

2a , $Xj + X2 > 0 . 


We can think of f(c^jc,t) as stp 

(-a(t) ,b) 

becomes that of Fig. 1.54, 


(£ £) • 


The system of Fig. 1.17 now 



where G(s) = Now f(«,t) is a positive operator. If we have a 


s ^+1 


resistance r ^ ^ in series with the inductor, then G(s) = 


3s-H 

s^+rs+1 


is a 


positive real function, 2 md so by the Positive Operator Theorem we will 
have the required stability. As yet we have not proven stability for more 
general cases. 


§1.9 Practical Considerations 

We now discuss briefly some of the practical aspects of implementing 
the control laws we have been considering. Fig. 1.55 shows a schematic for 
one possible implementation of the type of regulator we have been discussing. 
The two-position control switch is effected by means of the return-path 
diodes and D^, the power transistor Tj^, and the driver transistor T^. 

The two resistors marked are used to obtain a measure of half of the 
output voltage V^. This and half the desired output voltage of OE are fed 
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into the operational amplifier network which gives k(Bs+l) (V -as) , i.e. 

. R, 

k(0Vo+VQ-oiE) , where u * ^2^2 “ 25”* comparator gives a digital 

2 

(0 or 1) output depending on whether the output of the operational ampli- 
fier is positive or negative. This digital signal controls the driver 
transistors through the digital delay, the purpose of which is to limit 
the switching frequency of the transistors, so that the overall efficiency 
remains high. The parameter 6 should be chosen to give a good transient 
response to changes in load. The easy extension of this method of imple- 
mentation to the fourth-order case should be clear. Note that when a load 
resistance is present we are more concerned with V. than with I , so there 
is no need to measure inductor current. 

In practice the operational amplifier will involve a lag term in the 

measure of (V^-OE) , so that the voltage at its output is more accurately 
l4*8s 

given by k (V q-oe) , where 6 « $. Alternatively, it may be necessary 

1 . . 
to introduce this lag intentionally by placing a capacitor m 

parallel with R^, where = 6, so that the operational amplifier network 

does not become too receptive to hign frequency noise. 

We must therefore re-examine the system depicted in Fig. 1.17, and 

ascertain whether this is stab.e when the forward-path transfer function 

is generalized to be of the form G(s) = — . A short calculation 

(s^+1) (6s+l) 

shows that the associated linear feedback system of Fig. 1.19 is stable 
for all k > 0 provided 6 < 6. We therefore know that a positive real 
multiplier Z(s) exists such that G(s) 2(s) is positive real. The first 
multiplier to try is the simplest, i.e., Z(s) = (1+ys) , as used in §1.4 (d) 
and §1.6. We find 
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G(s) Z(s) - . <PfA) 

(sVl) (6s+l) 


Re G(ju) Z(ja>) = 


1 - (By-6B»<5y)(o^ 
(1-u^) (1+6 V) 


f 


so Re G(jw) Z(jw) > 0 for all w if and only if BySS-fiy «= 1, i.e. 

1 ft 

Y “ g ' j ' g " " ; and Y ^ 0 if and only if 6 < 6. We thus find that the same 

type of multiplier as considered for the simple second order case can be 

used here to prove stability. The reason for this is, essentially, that 

in that case we had Re G(ju) Z(ju) » 1 for allu , i.e. we had some 

"room to spare". The important conclusion here is that since the lag 

term ((6s+l) in the denominator of G(s)) will always be present in practice, 

some form of phase advance ( (6s+l) in the numerator of G(s)) is necessary’ 

for stability, with B > 6. Indeed, an analog computer simulation shows 

that if B is reduced below 6, an oscillation will occur, at a frequency 

of the order of — ^ — . 

2ir/[c 

References (25] and [14] are recent piiblications summarizing state- 
of-the-art techniques used in practice in the design of solid-state power 
supplies. Both give actual design examples of second-order switching 
voltage regulators which are similar to that considered here? ([25], 

Vol. 2, p. 165; [14], p. 196). The circuit given in [25] is of the type 
shewn in Fig. 1.56, which will be seen to be very similar to that shown in 
Fig. 1.55. Frequency limiting is effected by the hysteresis in the Schmitt 
trigger, the circuit of which is shown in detail in Fig. 1.57. The Schmitt 
trigger can be approximated (very roughly, but adequately for our purposes) 
as being equivalent to a linear system with transfer function (6s+l) , 
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followed by a hysteresis device with two-state output. The important 
point is that since our analysis has shewn that $ must be nonzero > we 
see that the 0.15 yF "commutating" capacitor in the Schmitt trigger circuit 
plays a crucial role in determining stability and ensuring a reasonable 
transient response, since it is responsible for the phase advance term 
(0S+1) . Both authors (25 ] , (14 ] give considerable attention to other less 
important design features, but do not even mention stability considerations, 
let alone explain the presence of this capacitor. Furthermore, neither 
author considers the possibility of using a fourth- or higher-order filter, 
which we have shown to be just as easy to implement and steibilize, while 
having superior design advantages. 

Reference (15 J also considers a second-order switching regulator, in 
which the Schmitt trigger is considered to have no dynamics, i.e. it 
cannot xntroduce a stabilizing phase advance factor (Bs+l) . The regulators 
considered in (15] are assumed to have a small resistance in series with 
the output filter capacitor. This will introduce the desired phase advance 
term in the forward-path transfer function. We might call this an 
"accidental" stabilization, which is another example of dissipation aiding 
stabilization. Clearly there is no control of the transient response in 
this design. 

In conclusion we note that all of the various methods of stability 
analysis which we have considered in Chapter 1 are needed for a lull 
understanding of systems of the type we ha^'e been considering. 
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CHAPTER 2 

POSITIVE OPERATORS AND DISSIPATIVE SYSTEMS 
§ 2.1 Introductio n 

The intent of this chapter is to outline the relevant background 
material for Chapter 1. In § 2.2 we follow the development of 
reference [ 34 ] in providing a simple proof of the Positive Operator 
Theorem, which we believe is a useful theorem, though so far has been 
applied relatively little. Fox more extensive and rigorous treatments, 
the reader is referred to [ 34 ] and ( 42 ] . In § 2.3 we discuss concepts 
of Dissipative System Theory following the development in the recent two- 
part paper by Willems, [ 35 ] , [ 36 ] . In § 2#4 we address the problem of 
obtaining state-space realizations for transfer functions of the O'Shea 
type, as introduced in § 1.5(d) and (e) , Theorems 1.4 and 1.5. 

§ 2.2 Positive Operators 

(a) Operators, and functions of time 

We consider functions of time on the interval 0 ^ t < The 
functions will all be real-valued, though our statements and theorems are 
easily generalized to the case of vector-valued functions of time. An 
operato r F maps a function of time x(t) into another function of time y(t) ; 
we write y ® Fx. Usually we think of these two functions x and y as the 
input and the output , and of the operator as an input-output sy r*"em. An 

-89- 

PKJdCiliJlNG PAGE BLANK NOT FIIiQlX 


I 



I 


-90- 


operator nay be specified by a characteristic graph, a convolution integral, 
a transfer function, or other means. We assume that all operators con*> 
sidered are causal, i.e. present and past values of the outputs do not 
depend on future values of the inputs. ^2 operators and a 

is a real number then the operators and aP^^ are defined by 

(Fj^ + F^)x “ Pj^x + PgX and (aF^)x = a(Fj^x) respectively. The operator 
PjFj^ is called the composition of P^ with P^ and is defined by (F^F^^jx = 
F^iT^x). In general F^F^^ ^ F^F^. If F 2 F^^ ■ F^F^ then we say that F^ and 
F^ commute . An operator F is linear if F(Oj^x^ + ^ 2 * 2 ^ ~ ^l^^'l * 
holds for all Xj^, x^. The identity operator I is defined by lx » x, and 
the zero operator 0 is defined by Ox = 0. An operator F is invertible if 
there is another operator F ^ such that FF ^ » F ^F = 1. (^ 2 ^ 1 ^ ^ * 

An operator of importance to us is the truncation operator P^, 
(more usually called a projection operator) , defined by 


(P^x) (t) - 


x(t) for 0 t ^ T < 
0 otherwise . 


An operator F is causal if P^F commutes with P^ for all T, i.e. 

2 

for all T. (Note that P^ = P^) . This is equivalent to requiring that 

Px, = P„x^ P„Fx, * P„Fx- for all T, (provided PO = 0) . 

Tl T2 T1 T2 

We assume that the reader is familiar with the concept of a vector 
space . A vector space V is called a normed vector space if a map, called 
the norm and denoted by 1 | 1 1 , from V into the real numbers IR is defined 
on V, such that for all x, y C V '\nd a C R, 
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(1) llx|^ 0, “ 0 if and only i'.' x = 0 

(2) I |ax| I = |a| I |x| I 

(3) ||x+y|| $ ||x]| + ||y|j (the triangle inequality ) . 

A real inner product space is a vector space V with a map, denoted by 
<»> and called the inner product , from V x V into IR such tliat <‘or all 
X, y, 2 e V and a, $ e IR 

(1) <x,y> = <y,x> 

(2) <ax+by, 2 > <= a<x,z> + 3<y,z> 

(3) <x,x> ^ 0, = 0 if and only if x = 0. 

An inn»r product space is b normed vector space, with ||xl| = /<x,x> . 

The Cauchy-Schwartz inequality states that 

!<x,y>| < IWI . ||y|| . 

The inner-product space which we consider is the space of real-valued 
square-integrable functions of time, L^IO,**]. A function x(t) is in if 
J* X dt < For X, y e the inner product is defined by 

<x,y> ^ f xy dt, so that l|x|| = x^ dt* 

0 ’ 

There are many functions of interest to us which are not in for 

example the constant functions, or functions such as t and e^. In order 

to be able to handle these we introduce the extended space which 

consists of all those functions x(t) for which P^x e for all finit"* T. 

T 2 

This includes all furctiens y(t) for which t < ® -=>ly(t) | < » . 
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We assume henceforth that all operators considered satisfy FO « 0. 
The operator F is said to be bounded if 


Fx 


sup ' ' ' < » . 

11x11 

This supremum will be called the gain of F, and denoted IIfII. It is easy 
to show that ^ HFjII. \\f^\\, that Hp^xM $ HxH for all T, 

and that 1 Ip^fI 1 $ 1 1f 1 1 for all T. F is said to be Lipschitz continuous 
if 


x,yeL, 

d 

x*y 


.„r I |F«-Fy| I , . . 
I l*-y| I 


This supremum will be called the Lipschitz constant of F, and denoted 
1 |f 11^; it satisfies the inequalities just given for IIfII. 


(b) Positive Operators 

Let FtL^ L^f i.e. F maps L^ into L^. Then F is said to be positive 
on L^ if for all x e Ljr 

<x,Fx> » 0 , 

i.e. I xy dt ^ 0 where y ® Fx. 

•'o 

c.- ;2. 

1 


i 



^ I « 1 , . < 1 , 

I 

I 

fi' 
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I 

I 

* 

F is said to be strictly positive on if P - r)I is positive on for 
some T) > 0, i.e. <x,Fx> ^ ri||xlP. 

Suppose now that ^ 2 e' ^ said to be positive on 

if for all x C and all T < », <P^x,P^Fx> ^ 0, i.e. j xy dt ^ 0 

for all T, where y = Fx. F is strictly positive on if F - r)I is 

positive on for some ri > 0 . The relationship between positivity of 

an operator on and L. is simple: 

2 2e 

Lemma If F:L_ "*■ L. and P:L- -► L_, then F is positive on L. if and 
' ' '■ 2e ^e 2 2 2e 

only if it is positive on L^. 

Proof Since F is causal, for all x e L. we have 
' 2e 

<P^x,P/x> * <P,pX,P^FP^x> = <P^x,FP^x> , 

which shows that positivity on implies positivity on L_ . Now assume 

that F is positive on L^^, but that for some x C L^, <x,Fx> < 0. Since 

lim P X = X, this implies by continuity of <,> that for some T, 

’£■*<0 ^ 

<P„x,Fx> » <P^x,P^Fx> < 0, which yields the contradiction. Q.E.D. 

Note that the lemma still holds if "positive" is replaced by "strictly 
positive" . 

If Fj^ and F^ are positive and and are nonnegative real numbers, 
then *^2^2 positive. If F ^ exists and F is positive then F ^ 

is also positive, (since if y = Fx then <F ^y,y> * <x,Fx> ^ 0) . Note that 
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/ •T 

^ xy dt ^ 0 can be viewed as a statement about the 
correlation between the functions of time x and y, a point of view which 
is utilized in the proof of the O'Shea theorem (Theorem 1.4 ) , [ 12 1 • 

Examples of positive operators are (i) A positive linear gain, 

(Fx) (t) = kx(t) where k i 0, (ii) A first- and third-quadrant nonlinearity, 
(Fx) (t) * f(x(t)) where fff ( (T) ^ 0, e.g. the functions stp C, sgn (T, 
sat k O, defined in Chapter 1, (iii) A time-varying nonlinearity, (Fx) (t)= 
f(x(t),t) where at ( ff,t) > 0, as for example in § 1.8. All of these 
examples are memoryless operators? that is, present values of the outputs 
depend only on present* values of the inputs. An operator which is not 
memoryless is said to be dynamic . An important class of dynamic operators 

is the convolution operators, i.e. those for which Fx * y where y(t) = 

/ t 

q(t-T)x(T)dT. Usually these are defined by a rational transfer 

0 

function G(s) , which is a function of the complex frequency variable s, 

representing a differential equation relating y(t) and x(t) , and given by 

e g(t)dt. Such an operator will be positive on L_ if and 
0 2e 

only if G(s) is a positive real function , as discussed at the beginning 
of § 1.4(d). This requires that Re G(joj) ^ 0 for all u, and strict 
positivity requires that Re G(jw) ^ r) > 0 for all co, for some 
T) > 0. 

To indicate why this is true we make use of Parseval's equality of 

Fourier transform theory, ((34] section 1.3, [ 29l)» We extend our time 

interval of definition to (-",») by letting x(t) = 0 for -®<t<0. Now if 

x(t) e L-(-®,'”) let X(jw) = lim x(t)e ^*^^dt. Then X(jo3) e L* (-",") , 

^ •J>4«D -T * 
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and X(ju) is called the limit-in-the-mean Fourier transform of x(t) . We 


have the inverse transform x(t) * -t— lira J X(jw)e du). Parseval': 

2ir vj^ <-w 

1 /*co 

equality states that J u(t)y(t)dt « rr J 0(jw) Y(jw)da) , where 

overbar denotes complex conjugation. 

Suppose now that y(t) = j ^ 

•'o 

g(t“T)u(T)dx, so that 

Y(jO)) » G(jU)U(jW). 


Then J" u(t)y(t)dt » 

J* U(jw) Y(ja))do) by Parseval 

*^00 

It 

f U(jw)G(ja))U(ju)dui 
—00 

II 

|u(jto) l^G(ju)da) 

VI 

II 

f |u(jw) 1 ^ReG( jw)dii) since u(t) and 

-00 


y(t) are real. 


From this it follows that uydt ^ 0 for all u e L^tO,*) if and only if 

•■00 

Re G(jb)) > 0 for all u). 

Finally, another important class of positive operators are those 
obtained by the composition of a monotone nonlinearity F with a transfer 
function G(s) . For a specific nonlinearity F, such as f{0‘)s c® , it 
is not known at present how to determine conditions on G(s) for FG to be a 
positive operator. However, if we require FG to be positive for any 
monotone nonlinearity f with f(0) « 0, then the theorem of O'Shea (our 
Theorem 1.4) provides the answer. 
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(c) The Positive Operator Theorem 

Consider the feedback system of Fig. 2.1. 



Fig. 2.1 


The governing equations for this system are 


/ 


ui - y 

U2 + y 


2 2 


2 

1 


We assume that the operators and satisfy G^ 0 = 0, (i = 1, 2) . We 

call u^, u^ the inputs ; the errors ; and y^, the outputs . The 

inputs may represent driving functions, driving noise, or initial condition 

responses. Assume that for any u. , u* e L_ , solutions e, , e_, y, , y_ in 

1 2e 12 12 

L exist for the above equations, and depend on u. , u in a causal way; 

2e 1 2 

this assumption is called vell-posedness . 


We would like to know when the feedback system is stable in the 


sense that bounded inputs yield bounded errors and outputs. 
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I 
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Definition The feedback system of Fig. 2.1 is finite-gain stable if any 
inputs Uj^, e Lj yield e^^, Oj, y^^* yj e Lj' there exist 
constants < » such that for any u^^, u^ e 

llejl, lUjII, lly^ll, llvjll S pJIuilU . 


Theorem 2.1 ( t 34 ] , Chapter 4) . Suppose that is bounded and either 

Gj^ is bounded or IIp^G^x]! ^ T]||P^xIl for all x e for some r) > 0. 

Then the system of Fig. 2.1 is finite-gain stable if and only if 
||(I + G^Gj^)'^!! < « on Lj. 

The theorem holds if the roles of G. and G are reversed. 


Proof Since ~ \ ~ G^e^ ® ' S^°l®l‘'’“2^ ' ^ * ^2^1®1 

^2^^1®1^'*2^ ' ^ exists on and is causal, by well-posedness. 

Let Uj^, u^ e l >2 be given. Then 




NVi * V2V1 - V2'°iW II 



lIVlII * IIWtVi * ‘’tVt'Vi^“2>II 


$ 

IIVill * IIV2lla IIVi'i ■ ”t'Vi''“2’II 


$ 

ll«ill * ll«jll • INjIi • 

Now 

u 

HI 


So 

Pm®! “ 
T 1 

?j(I*G2Gj)-^I4G2Gj)ej 


B 

P^(I*GjGj)“*rj(I.GjGj)ej since (ItG G l"'- is causel 


I - I rt . .... I - 4 ^ , 1- > V *<*.■>'« .. (f 

I 


i 


J 


J 



1 


1 


1 


t 


I 


* 


] 


Thus ||P^e^l| 

$ 

Therefore, ||p^ej^|| $ 1 | (I+G2G^)"^| | ,| |u^| | + | | (I+G^G^) ‘^ | M | G^ | | . | | | | , 

for all T. Thus e L^, and the finite-gain condition is 

satisfied. Since ^2 satisfies these conditions. Now if 

I |g I I < » , we see that y. . e^ also satisfy these conditions. Alternative- 
ly, if hMp^e^M ^ l|Pjy2ll then 62EL2 with ||e2ll< n"^||y2ll. and y^GL2 also, 
since necessity of | | (I+G^G^) ^|| < for finite-gain 

stability follows from 

e^ = (I+G2G^) , if we take = 0 . QED 

Theorem 2 . 2 . Positive Operator Theorem ([ 42 ], [ 36 ] Chapter 4 ). The system 
of Fig. 2.1 is finite-gain stable if G^ and G^ are positive on and 

one of them is strictly positive and Lipschitz continuous on L^. 

Proof * Suppose G^ is strictly positive and Lipschitz continuous. We shall 
show that the conditions of Theorem 2.1 are met. For any x e we have 

2 

n||P^x|| $ <P^G2 X,P^x> by strict positivity of G^ 

$ ||PtG2x||. I |p^x| I by the Cauchy-Srhwartz 

inequality. 

Thus, I jp^G^xI I ^ ri| |p^x| I for all T. We therefore only need to show that 
I I I+G2Gj^)~^| I < 00 on L^. So let x e and suppose (I+G2G^)y = x. Then 
y c by assumption, and P^y + P^G2Gj^y = 
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||Pt(I+G2Gi)"^ P^(I+G2G^)eJ| 

I|(I+G2‘^i)”^I |. I I ' 
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Now < p^y + ^ P^G^y,P^y > + < ^ 

» 0 + n||p^,Gj^yII^ 

since Gj^ is positive and G^ is strictly 
positive. 

Thus n| |p,j.Gj^y| 1^ ^ ^ ^ ll^T^iyil 

Cauchy-Schwartz inequality, 
so llp^G^yll $ n-^llp^xll $ n’^llxl I . 

Therefore, G^y e and | jc^yl | < n“^|!xl|. Furthenaore, 

l|G 2 G^y|( $ liG^IMlG^ylN h'^llG^ll. I|x||. Now y = -G^G^y + X, 
so y e and | |y| | = I l-G^G^y + x] [ 

5 llG2GxyM + 

^ h"'llG 2 l|.I|x|| + ||x|l . 

Thus, since y « (l+G^Gj^)”^ x, we have j | (I+G 2 Gj^) ”^1 1 ^ n”^llG 2 l| + 1 < » , 
and the conditions of Theorem 2.1 are satisfied. QED. 

It is interesting to note that if in addition G^”^ exists, then 
G (I+G G ) ^ turns out to be itself strictly positive and bounded on L- 

X A X » 

([341, p. 39) . Since y^^ = Gj^(I+G.,Gj^)”^Uj^ when U 2 = 0, this says that the 
closed-loop system itself is a positive operator. We can interpret this 
in terms of passive electrical networks, [42]: Let Gj^ be the driving-point 

impedance of a passive two-terminal network, and let G 2 be a passive driving- 
point admittance. Then 

the driving-point impedance of the series connection, as shown in Fig. 2.2, 
and this is passive if G^ and G^ are passive. 
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Fig. 2.2 

It is also interesting to note that the requirements of Theorem 2.2 can be 
altered to strict positivity of and G^, but no boundedness. 


§ 2.3 Dissipative Systems 

In Chapter 1,§ 1.4(e),§ 1.5(e),§ 1.6, and § 1.7(c) we introduced and 
made use of the idea of a dissipative system. Dissipative systems are of 
interest in engineering and physics; typical examples are passive electric- 
al networks in which part of the energy is dissipated in resistors as heat, 
viscoelastic systems v/here energy is lost through viscoelastic friction, and 
thermodynamic systems for which the second law of thermodynamics postulates 
a form of dissipation leading to an increase in entropy. We use the term 
dissipative as a generalization of the concept of passivity, and the term 
storage function as a generalization of the concept of stored energy or 
entropy . 
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Let G be an operator » or eystem« with input u(t) and output ^(t ) , 
0 $ t < *» defined by the equations 


x(t) • 

I - a(x) 


x(0) 


“0 


The vector >^(t) is called the state vector for G; the space X to which 
x(t) belongs is called the state space for G. Let w(u»^ be a real-valued 
function of u(t) and y{t). Then G is said to be dissipative vith respect 
to the supply-rate w(u,y) if there exists a nonnegative function V(jO with 
V(0) ■ 0 such that 

^ - w 0 . 

V is called the storage function for G. The inequality V - w $ 0 is called 
the dissipation inequality , and it is easy to show (using for instance 
Theorems 6.10 and 6.15 of reference [27] that this is equivalent to the 
inequality 


/ w(t) dt ^ 

^0 

where *^^0* " — 0 ' -^^ 1 * “ Sj. 

Mow we introduce a quantity called the available storage ; it is a 
generalization of the concept of "available energy" or "recoverable work". 
For the network of Fig. 1.34 it represents the maximum possible energy 
available at the terminals » starting from some given initial condition. The 
available storage of the system G in state ^ is defined by 

h 

V^(Xq) ■ sup - / w(t) dt 
x(0)«x Q 


t 


1 



i 


I 


1 


1 


1 
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The supremum is taken over all motions starting in state x ^ at time 0, 
and all possible ^(t) . The available storage is an important function in 
determining whether or not a system is dissipative, as is shown by the 
following theorem: 

Theorem 2.3. [35] The available storage V is finite for all if and 
— ■■ ■ ■ a —0 

only if G is dissipative. Moreover, 0 ^ V t V for iiny storage function 

a 

V, and V itself satisfies the dissipation inequality, 
a 

The reader is referred to [35] for a proof. 

The state space X of the system G is said to be reachable from 
x^C X if for any x^ e X there is an input function u^(t) which will 
transfer the state of G from ^^to sc^. We shall assume henceforth that 
all storage functions V have the property V{£) ■ 0, i.e. £ is the point of 
minimum storage for G in X. 

Next we introduce another quantity, called the required supply; 
this is the minimum amount of supply which must be delivered to the system 
in order to transfer it from its state of minimum storage (the zero st^ltc) 
to some other given state. The required supply for the state Xj^ of the 
system G is defined by 

V^(x^) - inf 

V® 0 

x(0)-0 

The infimum is taken over all possible motions starting in state >^(0) ^ 0 
and terminating at time t^^ in state x^. We now have: 


1 

/ 


uy dt 


i 
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Theorom 2.4 [35] Let G be dif sipative with storage function V for which 

V(£) ■ 0. Then 0 $ V < V < V . Moreover, if the state space X of G is j 

** A IT 

If 

reachable from £ then < « and is itself a possible storage | 

■i 

function. 

For a proof of this theorem the reader is again referred to [35]; we show 
that satisfies the dissipation inequality in the proof of Theorem 2.S. 

From Theorem 2.4 we see that the storage fvinction V of a dissipative 
system always satisfies the inequality ^ V ^ , i.e. a dissipative 

system can only supply to the outside a fraction of what it has stored and 
can store only a fraction of what has been supplied. and themselves 

always satisfy the dissipation inequality, and hence are storage functions. 

However, not every function V^(x) which satisfies ^ ^ will be a 

storage function. It appears to be difficult to state other general proper- 
ties of the set of possible storage functions, e^'cept for its convexity; 

If and Vj are storage functions, then sr s + (l-alV^ for any 

0 ^ a $ 1. This follows immediately from ' issipation inequality. In 
particular, if the state space is reachable from 0 then aV + (l-a)V is a 

•“A JT 

storage function for any 0 $ a $ 1. 

As a consequence of the normalization V(£) 0 we obtain the follov.’ing 

expected relationship between positive operators and dissipative systems: 

Theorem 2.5 Let an operator G with input u(t) and output y(t) be defined 

X ■ f (x,u,t) 
y - g(x) 

0 € t < • 



1 
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«nd assume that the state space o£ G is reachable from 0. Then when 
jc Q • £ G is a positive operator if and only if G is dissipative with 
respect to the suppiv rate uy with a storage function V(£) ^ 0 for which 
V(0) « 0. 


Proof . If G is dissipative then when x. ■ £ we have 


J uy dt ^ y V dt « V(x(T)) - V(x(0)) 


V(x(T)) 


^ 0 , 


and thus G is a positive operator. Suppose now that G is a positive 


operator whenever x^ ■ O.* Then 


V(x ) - inf f 
0 

x(0)«0 

xTtj^)«=Xj^ 


/• 

f uy dt which is the 


required supply. 


Then V(xj^) i 0 since ^ uy dt 0, and V(£) « 0 by taking = 0. 

To show that V-uy $ 0 we observe that 


V(Xo) + / uy dt J V(Xj^) 

t t ^ 

inf / uy dt + / uy dt ^ inf / 

S*” 0 t 0 


X (0) “£ 


X (t 1 *x ® 
— 0 -0 

x(tj^)«Xj^ 


x(0)-0 
X (t^) «=x^ 


uy dt 


I 
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Fig. 2.3 

Pig. 2.3 illustrates this last Inequality, which we see is precisely the 
dissipation inequality. QED* 

Note that we could also make use of the available storage to define V, 
rather than the required supply. The advantage of the positive operator 
concept is that it is an input-output concept which does not involve 
introduction of state space notions. However, in this approach we are 
restricted to considering operators which map 0 into 0, which though not 
a serious restriction, does in fact mean that we know something about the 
internal properties of such an operator, i.e. it must start from a state of 
mini.'"jm internal storage; the dissipative property of a system is indepen- 
dent of its initial condition. In applications both viewpoints are 
essential. For instance, in Chapter 1 we saw that for the second-order 
systems e dissipative system characterization of operators was superior 
since it led to Lyapunov functions, while for the fourth-order systems we 
had to fall back on the positive-operator methods. 

If we write d ® V - w then d is called the dissipation rate . When 
d E 0 the system G is called lossless . An example would be the driving- 
point impedance of a network containing only inductors and capacitors. 
Energy storage and retrieval is 100% efficient for such a system. In 
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this case the dissipation Inequality becomes an equality, and the storage 
function is defined uniquely by »• V ■ ■ / w dt. 

It may happen that for a system which is not lossless; such a 

system is called quasi-lossless since it can be transferred between states 
with arbitrarily small dissipation if the input is suitably chosen. All 
first-order systems, (such as that of Fig. 1.32) are quasi-lossless. As an 
example consider the system of Fig. 2.4. where f(cT) is a function with 
f(0) =0 and 0f(0) >0. Let F(x) = f f(0)dO. 


U* X + X 


t(<r) f(x)=y_ 


Fig. 2.4 


If w » uy = (x+x)f(x) = xf (x) + £(x)x we have 


V (X ) - sup 
t^^O 

x(0)=x. 


-/[ 


xf(x) + Jdt 


[-/ 


xf(x)dt - F(x(tj^)) + F(x(0))J 


x(0)=x- 









(• 



-107- 


Now we can always choose u so that xf(x) is bounded. Thus* letting x(t^) 
“ 0 and tj^ 0 we obtain (x^) ■ F (x^) . Similarly 


V (x.) “ inf 
' ‘ tjSO 

x(0)=0 

x(tj^)-Xj 


I xf(x)dt + F(Xj^) 


» F(Xj^) . 


Let us consider an electrical interpretation of this when f(o) ■ a. In 
Fig. 2.5 u is the applied input voltage to a resistor-inductor network, 
and X is the resulting current, given by x + x ® u. If x(0) = 0 then the 
energy supplied to the network in the interval 0 < t $ t^ is 


1 

r 2 1 2 

ux dt = J x dt + — Xj^ , where Xj^ = x(tj) . 

° 0 ,2 

The energy stored in the inductor at time t, is E. = — x (t, ) . In order to 

12 2 1 

demonstrate the quasi-lossless property of this network we must exhibit a 

t 2 

u(t) for 0 $ t $ tj^ which will make the difference | ^ 


arbitrarily small. One such u(t) is the constant function u(t) 


(^1 
r 


x^d - e’^) 

for 0 $ t $ t , which gives x(t) = ' ■ $ x 

1 - 


1 ' 


2 2 

so that J X dt < tj^Xj^ , which approaches zero as tj^ 0. Thus, the 
0 


ia 



1H 


Fig. 2.5 
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most e££i.cient driving £unction for this network is an impulse, i<e. for 
high efficiency the energy should be delivered to the network in as short 
a time as possible. 

Consider now the linear system G described by 
X = A X + B u 

Z. “ X + ^ \i 

where x(t) e IR*' ; u(t) , y(t)eiR*'', and are constant matrices of 

appropriate dimensions. Assume that (A>B,C,D^) is minimal. We wish to 
know when this system is dissipative with respect to the supply rate 
w * u'^. If D + D' is invertible the evaluation of the available storage 
V and the required supply V reduces to an optimal control problem which 

1C 

may be solved by considering appropriate solutions of the matrix equation 

H A ®“£'> <£ + = 0. 



It can be shown [36] that this has a real symmetric positive definite 

solution if and only if G is dissipative. Then V (x) = ^ x'K x and 
1 + + - 

V^(jc^) = — x*K x,where 1C and K are solutions of the above matrix equation. 

If D + D' is not invertible then and K** are given by the limits as 

n 0 of and K ^ , which are solutions of the above matrix equation 
— h — n 

obtained by replacing £ by (D^ + q^) . We also have that G is dissipative 

with respect to w = if and only if there is a real matrix 2. “ 2.' ^2 







^ ^ 1 ’ ] 1 1 
!• i i ! . : . i 
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Moreoveri ^ x'2. 2 . * storage function if and only if 2. satisfies this 

inequality. In particular and K satisfy it, and every solution £. 
satisfies 0 < K~ $ $ K^. 

§ 2.4 Realizations for O'Shea Functions 

Suppose we have an operator F with input u(t) and output y(t) defined 
on 0 $ t < « by the equations 

2 . “ ^ 2 . ^ 

z ® £ 2. 

y = f(z) 

x(0) = 0 

where f(0) is a monotone function with f(0) >= 0, as depicted in Fig. 2.6. 

As stated in Theorem 1.4, we know that F is a positive operator if and only 

00 

if Z(s) = [c(Is-A)’’^b]”^ = g„ + Ys - g(s), where Y ^ 0, g(s) = / g(t)e"®^dt, 
g(t) > 0, and g^ ^ g(0) = j g(t)dt. If we also know that f( (T) is odd, 
i.e. f(0l « -f(-cT) t then we can relax the conditions on 2(s) to require 




Fig. 2.6 
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[32] « [12], and we call a transfer function for which z(s) neets 

these conditions am 0*Shea function . 

Given am O'Shea function z ^(s) we would like to be able to construct 
a realization (A,b,£) of z(s) , as in Fig. 2.6, which will allow us to 
write down a storage function for F as a function of the state vector 3c. 

As yet a general solution is not known to this problem, however, we have 
the following theorem, (cf. Theorem 1.5) : 

Theorem 2.6 Suppose that -A(t) is a hyperdominamt matrix for all t, 

^' = [0 0 %«.*0 Y] , Y > 0 , £ = [0 0 0 1], and f is any monotone non- 

linearity with f '0) = 0. Then the operator mapping u(t) into y(t) defined 
by 

X = ^ 21 

z “ 21 

y = f(z) 

is dissipative with respect to the supply rate uy; with storage function 
1 V rZ 

^^^21^ “ T \ F(x.) , where F(z) = j f(o) do. Furthermore, if f is odd, 

i=l ^ 0 

then -A(t) need only be dominant. 

A matrix M whose ijth element is m. . is said to be dominant if 
n n ^ 

w. . ^ I 1^. . 1 ni, . ^ y |m. .1 for all i,j; that is, for each row and 

3«1 i«l 

column the on-diagonal elements are larger than the s\m of the moduli of 

the off-diagonal elements. M is said to be hyperdominant if it is dominant 

and all the off-diagonal elements are negative; i.e. M is hyperdominant if 

when i 2, and I m. , ^ 0 and ^ m . ^ 0 for all i,j. 
i=l j=l 
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Proof of Theorem 2 . 6 j 


*‘i’ ■ X Jj S*i - J i!, i ““'f 


where 


tl(x-) f(x.) ••• f(x )J 
X. i n 


ThuH V(x) = Y V A X + r V bu 

— A — X A — X — 

*-• V V Ax ^ uf (x ) 

A ~>i~— n 

“ ^ V ^ A X + uy 
Thcruforc: V-uy«^V^^^ . 

Now, it can be shown (( 33 ], ( 34 ] Chapter 3 ) that If f( 0 ) is any monotone 
function with f( 0 ) « 0 then is ncqatlve for all x if and only if 

is hypordominani. , and that if f( 0 ) is any odd monotone function thru 
nrqativo for all if and only if in dominanl . n\:u 


Ni'ces 5 ^ary and sufficient conditions on A for A x to be nccjalivt* foi 
all X, when f is a particular qiven function an* not known in cp^nnal. 
However there is one function f of inU^rest (pai t icularly with leqanl t< 
applications of Chapter 1 ) for which thc'se conditions arc known: this i.s 
the sgn function. I.et a row-dominant matrix M be one for wliich 


^ 2 ^ / nt. . for all i , 

A. ij 


sgn X, 


stj X 


and if X ^ x^ then let sgn x = sgn x and stp x stp x. 


sgn X 


stp X 


wc* have: 


i 


j 









-n?.- 

Throrom 2.7 x'M sgn x ? 0 for all x if and only if M is row-dominant, 
and (fj2£ x_) ' M )c 0 for all >t if and only if M is column-dominant. 

Proof . Sufficiency: x'M sgn x = *”12 *2 "'in*’*’” ^n^ 

+ *2 ^"*21 +••• ) 
4- • • • 

2 

= Im^j^Ccgn x^) + x^) (r.yn x^) x^) (s<jn J 


(assuming x^ fi 0) 


^ 0 if M is row-dominant 


Necessity: Suppose row 1 has m ~ i |m. .| < 0 , (the proof for the caiie 

11 ID 

•■“I l"'■•l ^ 0 being similar). Let x = (1, - T] sgn m. 5, *•* * “ *1 "'i.J * 

3L X . • XX X ~ X 1 1 

Df*! 


T\wn 5*.qn x -• (1^ -sgn m. -sgn m. -sgn m_ )*. 

‘ 1 ^ ■ X 3 


12 

n 


In 


Thus, x'M sgn x = m - T |m. . | - n sgn m (+m +••••• +m ) 

— * 11 13 * 1^ — 1^— 


j-2 


< 0 for rj sufficiently small - 


QKV. 


0n(‘ can also show in a similar way that a sufficient (but not nocosstuy) 

condition for x*M stp x to bo nonnegative for all x is that m. . ? P ^ 1^ .:! 

-■(-=,bl b a ' “ 

for each i, whore p is the larger of ~ and — . 

Lot us call a realization (Afb,£) of H(s) which meets the conditions 

of 'J’hcoiom 2.6 a dissipative (hyper ) dominant real iznt ion . Though it ir. 

not known what the necessary and sufficient conditions arc for a given H(r) 

to have a dissi)'ative (hyper) dominant realization, clearly it is necessary 

for H(s) to be an O'Shea function: H(s) = (g^ + Ys-g(s)) ^ with y ^ 0, 

I conjecture that this is also a sufficient condition; i.c. the conjecture 

in that every O'Shea function of the type (g + ^ with y ^ 0 

>/ 


\ 


k 
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has a time-invariant dissipative (hyper) dominant realization. Let us 
further discuss this question. 

Firstly, it is easy to show that the conjecture is true for a 
second-order O'Shea function, and that such a function need not have real 
poles. Secondly, it is always possible to find a realization (A,b,£)with 
b' = (0 ••.-O X] and £ •* [0***0 1] for a given transfer function H(s) pro- 


vided that Lim s H(s) “ X ^ 0. For, let 
s-> “ 


R 2 

^n-1 ° 

s" + Pq 


where q , ■ X 
n-1 


Then we have the st^ ..^rd controllable realization ([7] Chapter 17) 
H(s) = h(^s-F) where 




1 






A 
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Then we have the realization H(s) ■ £Us-Aj^) "^b where 


K-1 


- KL » b - L SL - 


0 

Vi 


c - h p“^ - [0*..0 11 


Now thus obtained will in general not be hyperdominant ; the conjecture 

is, however, that if H(s) is O'Shea then there exists an R such that 

»1 

® hyperdominant, R b ■ b , and c R = c . Note that if 


H(s) * (g^ + Ys-g(s))”^ withYy< 0 then Lim sH(s) « For R to satisfy 

s -► » 


.-1 


R b > b and c R 


-1 


c we must have 


^11 I a 

■rr 


In § 1.5(e) we considered an example of this approach for a third-order 
O'Shea function. 

Now let us consider another approach. An interesting subclass of the 

O'Shea functions is the class of RC impedance functions Z(s) which have 

Z(«o) ■ 0. A general RC impedance function Z(s) can always be expressed 

r ^i 

([2] , [16] , [42] ) as Z(s) * g„ + i — r — where g., s. ^ 0 for all i. If 

° i-1 ® i ^ ^ 

g^j = 0 then Z(«o) e o. An RC impedance has poles and zeros alternating on 
the negative real axis, and a Nyquist locus which lies inside a circle in 
the right half plane. It is proven in reference [42] that such an RC 
impedance is an O'Shea function. Now the four canonical RC network 
realizations of an RC impedance are the Cauer I, Cauer II, Foster I, and 
Foster II methods. Of these the Cauer I method is of interest to us 


i 


I 



here , because It leads to a tridiagonal realization which In some cases 
is hyperdominant. A tridiagonal realization (^,b,£)is one for which the 
matrix A is tridiagonal « i.e. the only nonzero entries are on the 
diagonal and inunediately above and below it: a^^ ■ 0 unless 

j e *^J.-l,i,i+l} . Consider the network of Fig. 2.7 w_th input current 
I and resulting terminal voltage V. 


We have 


C V ■ 
11 


V -V 
2 '^1 


V -V 
3 2 


V -V 
2 1 


C ,V , 
n-1 n-1 


V -V , V -V * 
n n-1 n-1 n-2 


- V -V , 

CV . ” . 1 

" " Vi 


A 

Letting x. * V. /c7 , u , y ® , we obtain 

’’ 1 i 1 n n 

I X “ h X + ^u 

I y ■ £ 2L 


Rn-1 Rn-2 



♦ 

r- ^ 

♦ 

r— WWW ^ 

♦ 

♦ j 

-.Vn - 



L - - ■ . -i 



Fig. 2.7 
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where 



b’ “ [0 0 ... 0 1], c * [0 0 ... 0 1], d * R C . Now -A is hyperdominant if 
f n n 

and only if « C for all i. Thus, if 2(s) is an RC impedance with Z(»)«0 

which can be realized by a Cauer I network with all capacitances equal , 

then Z(s) has a tridiagonal dissipative hyperdominant realization. However, 

a general RC impedance cannot be realized by a Cauer I network with all 

8+1 

capacitances equal, (For example — cannot). For n $ 3 it appears 

s'^+Ss+l 

that every RC impedance does indeed have a tridiagonal dissipative hyper- 
dominant realization; furthermore one can show that every such realization 
for n ^ 3 is an RC impedance. 

In obtaining a Cauer I network realization of a given Z(s), Z(s) is 

expressed as a continued fraction expansion. For instance, for 

2(s) « — we write 

SV3S+1 


i 





Fig. 2.6 

To obtain realisations for O'Shea functions in tridiagonal form we make 
use of analog computation symbols as defined in Fig. 2.9. (We leave 
Initial conditions unspecified) . 

y »-/(x^ + X2)dt 
y «-(x^+X2) 

y sOCX 





Fig. 2.9 
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NOW a transfer functior. G, (s) ■ — = 

1 as+f(s) 

realized by the configuration of Fig. 2.10. Similarly a transfer function 
" b+f(s) 


with input u and output y can be 


can be reil.lced by the configuration of Pig. 2.11, 




Now given any transfer function Z (s) with Lim sz(s) ft 0 we can always 

g ce 

obtain a tridiagonal realization (A,b,£) with b' » (0***0X1 and c « I0***01] 
by "“xpressing Z(s) as a continued fraction expansion and representing this 
expansion by a succession of connections as in Figs. 2.10 and 2.11. To 
illustrate this method consider the RC impedance 


s + 4s + 


n 

i 


s^ + 6s^ + 11s + 6 




1 , 


i 


1 


J 



Using the connections of Figs. 2.10 and 2.11 we now obtain the represen- 
tation of Fig. 2.12, in which we have made the three integrator outputs 
proportional to the state variables ^2 '^3* ^ labelling of inte- 

grator outputs gives rise to a tridiagonal realization. 



Fig. 2.12 


i 


1 



The conditions which a and @ must satisfy for this to be a hyperdominant 
realization are 




i 
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M 



This corresponds to the Cauer I network realization of Fig. 2.14. 


R.-I 


>,*re 


vf ; 


X Cg * 2 ✓ 

— — - — ^ 

-c .5oj 

3 2 


Fic 

[. 2.14 


In a realization of a given 2 (s) by this method, if we let = ^i^i 

is the output of integrator i and is some positive constant, then the 
A matrix obtained will have diagonal elements negative and off-diagonal ele- 
ments positive if and only if the coefficients in the continued fraction 
expansion are all positive. If this is so then 2 (s) can be realized as a 
(Cauer I) RC network jmpodance. 

Consider now the O'Shea function which we obtained in § 1.5(d): 
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G(s) 


. 6 ^ + 1.45s + 0.5 
+ 200s^ 4- 96s + 4 



This gives us the representation of Fig. 2.15, from which we obtain the 
state equations 



This realization cannot be made hyperdominant for any values of ot and 3; 
we need to label the integrator outputs of Fig. 2.15 with a general linear 
combination, i.e. we must label the outputs of the lower two integrators 
(otXj^ + Bxj) and (Yx^^ + 6 x 2 ) » then look for suitable .values of a,6,Y,6. 
This, however, is essentially the same task as finding a,b,c,d for the 
matrix R of § 1.5(e). 


I 


i 


i 
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CHAPTER 3 

SWITCHED ELECTRICAL NETWORKS AND BILINEAR EQUATIONS 


§ 3.1 Introduction 


A linear dynamical system can be described by a set of equations of 
the form 

^ ■ A X + B u 

( 1 ) 

Z. “ £ 2L 2. H. 


where u is the vector of inputs or controls, ^ is the vector of outputs, 
and 3 c is the state vector for the system. We saw in Chapter 1 that 
equations of this form describe the behavior of certain types of electri- 
cal networks: in Chapter 1 these were voltage-conversion networks 

operating from a battery with zero internal impedance, and we were adsle 
to exploit the fact that the resulting evolution equations were linear, 
for this led to a feedback system with a linear operator in the forward 
path. However, the class of systems describable by linear equations is a 
restricted one; most of the systems encountered in the field of electrical 
power processing (i.e. DC-DC conversion, DC-AC conversion, etc.) cannot be 
described by linear equations. In this Chapter we address ourselves to the 
equation "What kinds of state question arise in the description of power 
processing systems?" Having answered this question in § 3.2, we then ask 
in § 3.3 "What statements can ve make about classifying such systems: 
what canonical forms for the state equations do we have?" In § 3.4 we 
outline the role played by Lie groups and Lie algebras in characterizing 
these systems, and in § 3.5 we give network examples. 
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§ 3.2 Bilinear Equations 

Consider an electrical network composed of resistors, inductors, cap- 
acitors, transformers, batteries, current sources, and ideal switches. 

With these components, we can model the essential features of power con- 
version networks. The state of this network will be a vector j^(t) in IR*' 
whose instantaneous value represents the inductor currents and capacitor 
voltages; usually we take the state variables to be scalar multiples 
of these currents and voltages. 

If the topology of the network is fixed, that is, the switches are all 
held in one set of positions, the state evolution equations will take the 
linear form (1) given in § 3.1. The reason for this is that capacitor 
voltages and inductor currents obey linear first-order differential 
equations. The matrix A of § 3.1 will be a constant matrix v;hose 
eigenvalues represent the natural frequencies of the circuit. These natural 
frequencies are determined by the component values and the topology of the 
network, which will be changed if the switch positions are changed. 

Consequently, if '*'® have a network in which the switches are considered 
as controls, with the position of switch i being given by u^ where u^ = 0 

or 1, then the matrix A will be a function of the u^'s. The resulting 

state evolution equations take the form 

(2) s * “ u,A, + U-A- +**‘+u A )x + (b,u, + b.u_ +***+b u ). 

\ — “V 2“2 — “1 1 —2 2 — tn in 

As a simple example consider the regulator of § 1.7, shown in Fig. 3.1. 




^ l '^2 " ^1 


We have 




I 

I 

I 



Pig. 3.1 

Thus, letting and * V 2 (so that j J^'x is the stored 



i.e. X = (2^ + \iA )x + bu , 

When R » 0 then 2 state equation is linear; in terms of the 

network we see that when R ■ 0 the natural frequency is unchanged by the 
switch position, since the dyn£unic impedance of the loop including and 
is the same for both switch positions. This example is of the simplest 
kind, though of consider6d>le practical significance. In § 1.6(a) we 
briefly considered higher-order regulators of this type. 

Now any system which can be realized by a set of equations of the 
form (2) , i.e. 
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can be realized by a set of equations of the form 



To see this we let A ^ be obtained by adding an extra row and column to 



We are therefore interested in equations in the form of (4). This 

is called the bilinear form. This form is even more general than we have 

just shown it to be, for it is shown in reference [9] that if 

P(x) = P (x, ,x_ , • • • ,x ) is a multinomial expression in the varied^leg 
1 ^ n 

, then any system of the form 



can also be realized in the form of (4) . Thus nonlinear output maps can 
be reduced to linear forms provided that they are of the finite power 
series type. Again ^ this is done by extending the dimension of the state 
vector. For example^ the equations 

( X *= -X u 



can be written as 


1 
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This technique for handling nonlinear combinations of the state variables 
is of more than passing interest to us. Consider, for instance, the 
voltage regulator we studied in § 1.4. 

There we had 

X + X ■ u 

y ■ Bx + X 
u = -stp y . 

We foxmd that once the state reached the chattering region of the switch- 
ing line near the origin, the settling time was proportional to $. Thus, 
near the origin we want 0 small, while for large values of ||x|| we do not. 
In order to obtain an improved overall transient response we might 
therefore try a feedback strategy of the form y = x^ + x. We can now put 
the equations 

X + X - u 
• 3 

y “ X + X 
; u “ -Stp y 

in the bilinear form (4) . It must be pointed out that this technique will 
lead to state equations of large dimension; in this regulator example the 


state is 10-dimensional: 



-130- 


X. .[l 


XX 




XX 


k^l 


For a power conversion network with state evolution equations in 
the form of (2) or (4) , it will sometimes be the case that some of the 
variables u^ are restricted to being functions of the state. Such a 
situation may arise if there are diodes in the network, as for example in 
Fig. 3.2, for which 



Pig. 3.2 


we have 


where 


u 


jl . Xj < e 
lo , x^ ^ e 


In this thesis the only type of switch that we consider is the ideal two 
position switch depicted in Fig. 3.3. 



Fig. 3.3 




We use this switch In obtaining "well-posed” network models for power 
conversion systems; by a well-posed network we mean one for which capacitor 
voltages and Inductor currents are continuous functions of time. Fig. 3.4 
is an example of an ill-pcsed network (cf. Fig. 1.5). This is not well- 



Fig. 3.4 

posed because an infinite voltage would be developed across the switch if 
it were to be opened when the inductor current I was nonzero. By pro- 
viding an alternate path for the inductor current when the switch position 
changes, as in Fig. 1.5, the current will be a continuous function of time. 

Implementation of an ideal two-position switch will in general 
require the use of two transistors (or thyristors) and two "freewheeling" 
diodes, as for example in Fig. 3.5 which shows a scheme for implementing 


0 . 



Fig. 3.5 
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the regulator of Fig. 1.5 (or of Fig. 3.1 with R - 0) . When u - 1 tho 
transistor is turned *’on” (i.e. acts as a short circuit) while the 

transistor is turned "off” (I.e. acts as an open circuit) . In this 

condition^ the network of Fig. 3.5 is equivalent to that of Fig. 3.6, in 
which the inductor current 1 flows through when 1 > 0 and back through 


D| 



Fig. 3.6 

whenever I < 0. When u ■ C. is "off" and is "on", and then the 
network is equivalent to that of Fig. 3.7, in which the inductor current 
flows through D when I > 0 and through T when I < 0. 

» A 


T 



Fig. 3.7 


1 


1 


T 


T 


T 


T 




i 


-133- 

In practice « when conMtructlng a regulator of this type> it is 
usual to omit the transistor T (as in Figs. 1.55 and 1.56). This still 

m 

yields a well-posed network# however the analysis of Chapter 1 will 
apply only if the condition {u ■ 0# I < O) nev«ir arises# which will be 
the case if the load current is l*rge enough to ensure that at ail times 
I > 0. 

In conclusion# we see that the state equations for power conversion 
networks will be of the bilinear form (4)# with some of the restricted 

to being functions of the state# while others can be chosen freely. 

§ 3.3 Canonical Forms and Equivalent Systems 

Since a wide variety of complex electric power conversion networks 
have state equations of the form of (4) # we are particularly interested in 
classifying bilinear systems. We would like to be able to determine when 
two electrical networks which are topologically different have similar 
dynamical characteristics; not only would this be conceptually helpful# 
but a feedback law vAiich was found to be suitable for one network could 
be translated into a suitable feedback law for the other. Reference [9] 
describes some recent results which answer the question of when two bilinear 
systems are dynamically similar. We shall now briefly out ne these. 

Suppose we have two systems of bilinear equations: 

(-|>* Ji "A) - 

X. ■ C X 

H * 




I 


1 


Let u ■ (u, , u , ‘“u ) be the input vector, and y the output vector. We 
say that (4) and (6) are equivalent if there exists a nonsingular matrix 
P such that F ^ - P A ^ p”^ for 0 « i $ m and H *= C p"^. Our reason for 
this definition is that if we are (4) , and we let x.# then we 

obtain (6) with H “ 

are realizations of the same input-output map. We would like to know 
when the condition that (4) and (6) realize the same input-output map 
implies that they are equivalent. We shall answer this question below in 
Theorem 3.1, which is similar to the well-known result on the equivalence 
of realizations of a linear system ( [7] , section 18) . 

We call a realization in the form of (4) irreducible if there is no 
invertible matrix P such that for 0 i $ m 



where the are square matrices, all of the same dimension. That is, for 

no choice of basis is the realization in block triangular form. Otherwise 
it is called reducible . It is a fact (19], Theorem 3) that every bilinear 
realization (4) is equivalent to one in which the matrices A ^ are in block 
triangular form, with the diagonal blocks being irreducible, thus: 




I 


1 


1 


1 
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We shall call x ^ an equilibrixim state of the bilinear system (4) if 
0. This is the same as requiring that )Cq be an equilibrium 


solution of (4) when u^ = 0 for all i. 
We have 


Theorem 3.1 ( (9] , Theorem 8) : Suppose that we are given two realizations 

of the same input-output map 

m 


(4) 


x(t) = 


( 6 ) 


^0 " 

l 

i«l 

u,A . 
•-0.— 1 

C X 


x(0) 


m 



1 

i=l 

u.F. 

i-n 

H z 

# 

2(0) 


-0 


-0 


Let Xq and be equilibrium states. Suppose that both realizations are 

controllable in the sense that the set of states reachable from or 
z_Q is not confined to a proper linear subspace of the state space for 
X or Suppose also that both systems are observable in the sense that 
any two starting states (not necessarily and £^1 can be distinguished 
by means of the measurement of for a suitable choice of £. Then the two 
realizations are equivalent. 


§ 3.4 The Nature of Solutions for Bilinear Equations 

For the linear equations (1) given in § 3-1 it is well known that the 
explicit solution for x(t) in terms of u(t) is given, for all t, by 


i 


1 
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f 
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(7) I x(t) - e^Vo) + j u(T)dT . 

0 

This formula for x(t) is sometimes called the Variation of Constants 
Formula^ and is useful not only for explicit calculation of solutions, but 
also for determining properties of these solutions. 

It is not possible to write down an analogous explicit solution to 
equation (4) for all t. However, as shown in [8] we can study the intrinsic 
properties of the solutions of (4) by using the tools provided by the theory 
of matrix Lie Groups and Lie Algebras ([20], [28], [8], [3], [30], [10]). 

To this end, we now introduce these concepts. 

Let P. denote the set of real nxn matrices; it is a vector space of 
dimension n . A Lie algebra L in is a vector svibspace of which 

has the property that if X and Y belong to L, then so does [X,Y] = XY-YX. 

We call [X,Y] the commutator or Lie product of X and Y. As examp'' es of 
matrix Lie algebras we have: 

Example 1 : The set of all real nxn matrices IR*'**' is itself a Lie 

algebra, sometimes called the general linear Lie algebra and denoted gH{n,'R) . 
Example 2 t The set of all nxn real matrices with zero trace is called the 
special linear Lie algebra, and denoted sl(n,lR). 

Exeunple 3 : The set of all nxn real skewsymmetric matrices, i.e. those 
which satisfy X “ called the orthogonal Lie algebra and denoted 

o(n,E) . 

Example 4 ; The set of all 2n x 2n real symplectic matrices ,i .e. those which 
satisfy X'J + J X = 0 where 




, ,,-i 


.J. 
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I is called the symplectic Lie algebra and 


denoted sp(n,!R) . 

Example 5 ; Consider the three matrices 


R = 

—X 


0 0 0 
0 0-1 

0 10 


R 

-y 


0 0 1 

0 0 0 

-10 0 


R - 

— z 


0 -1 

1 0 

0 0 


These are a basis for oO^IR) , and we find that 


[R,R ]«R , [R,R ]*R , [K ]=R 

— y'— z' —X — 2— X — y — X— y — : 


Example 6 ; The affine algebra of the line, aff (1), consists of all real 

and 




1 

- 



matrices of the form 

u v 

. The two matrices X = 

1 C 




0 0 


o 

o 

Y = 

0 1 

form a basis, and 

we find that ® X.* 



O 

o 

1 






Given an arbitrary subset of we can add additional elements to ir 

so as to imbed it in a Lie algebra. To obtain the smallest Lie algebra 
which contains a given set N we first add to N all linear combinations 
of elements of N so as to obtain a vector space Then we form all pos- 
sible Lie products of elements in to obtain a set . We add this 

to and obtain . If N^*is not contained in we 

+ [N^# N^] , etc. This process must stop in a finite number of 

steps since otherwise at each stage we increase the dimension of the 

2 

vector space by at least one, and the dimension cannot exceed n . We call 
the Lie algebra thus obtained the Lie algebra generated by N , and we 


-J 


,,L «... 




I 
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denote it . For example « in Example 5 above the Lie algebra 

generated by any two of R ,R , R is o(3,IR). 

X y "" z 

Let L be a Lie algebra. The set [L,L] of all possible Lie products 
of elements in L is called the derived algebra of L and is denoted L*. 

The derived algebra of L' is denoted L". Continuing, we have the 
derived series ; 

LDL* DL-D L^“* D D ••• 

The Lie algebra L is said to be solvable if l^^ is zero for some m. We 
shall call L simple if L* » L. For example o{3) in Example 5 above is 
simple f while if L = aff(l) as in Example 6, then L* is spanned by Y and 
L“ is zero, so that L is solvable. It is a fact that any Lie algebra can 
be decomposed into the semidirect sum of simple and solvable parts, 

[28]. 

For each X in the Lie algebra L we define the operator ad by 

X 

ad Y » [X,Y] for all Y e L. 

Powers of ad are defined by 

adJJ Y = [ X,[^,y ♦,[X ,Y^] ...]] 
n tiroes 


Using this notation, we can state the following result, which we use in 
§ 3.5: 

Lemma (Baker-Hausdorff) : If X,Y are elements of the Lie algebra L then 

e^- Ye“- e L and Ae~- = Y + [X,Y^I + i [X, tX,Y^J ] 
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^ 
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For a proof 

of this 

the 

reader is referred to [19] . (to see how the scries 

is obtained 

let 




f(s) 


sjf -sX 
e — 

Then 

f(s) 

s 

sX„ „ -sX s^ „ -sX 

e— XYe “ 



8S 

— c V 

e -[X,y]e - 

Similarly 

f (s) 

S 

SX ,n „ -sX 

e - adj^ y e - 




“ 2 

Then 

f(s) 

s 

f(0) + sf • (O' + -|^ f"(0) + • •• 




2 



85 

Y + s[X,Y] +-|j [X, (X,Y)] + ••• 

and the result follows on putting s » 1.) 


Having introduced matrix Lie algebras we next introduce the con- 
cept of a matrix Lie Group. If M is a set of nonsingular matrices in 

1r”^, we let {m}_ denote the multiplicative matrix group generated by M, 

G 

i.e. it is the smallest set of matrices in R*^’**^ which contains M and 

which is closed under multiplication and inversion. If N is a linear 

subspace of then let P be the set of all matrices of the form 

Xi X, 2 d 

e e • * • e” where X . e N for each i and p = 0, 1, 2, ••* . P contains 

X. tr X. 

no singular matrices since for any matrix 2i^,det(e ) » e > Of 

([7], section 4). Since it is closed under multiplication and inversion, 

P is a group, and we write P * {exp n}_. It is an interesting and useful 

G 

fact that (exp n}_ = {exp {n}_}_ ([8] Theorem 1.) 

G AG 

If L is a Lie algebra then we call {exp l}^ the Lie group associated 
with L . For a full treatment of the relationship between Lie groups and 
Lie algebras the reader is referred to (30) and [10] . We now give the 
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Lie groups which are associated with the six Lie algebras given above: 
Example 1 ; The general linear group GL(n,!R) consists of all real nxn 
invertible matrices. 

Example 2 : The special linear group SL(n,IR) consists of all real nxn 

matrices with determinant 1. 

Example 3 ; The orthogonal group 0(n,lR) consists of the real nxn matrices 
which satisfy )C'X^ = The spe cial o rthogonal group SO(n) consists of 
all matrices in both 0(n) and SL(n,IR). 

Example 4 ; The symplectic group Sp(n,1R) consists of the real 2n x 2n 

matrices which satisfy X' J “ I, where 

0 -I. 

L 2, 

Example 5 ; The group S0(3) consists of all r^al 3 x 3 matrices which 

have determinant +1 and satisfy X^'J^ = 

Example 6 t The affine group of the line consists of all real 2x2 matrices 

[a b ' 

whose second row is TO 1] . Any such matrix is of the form 

0 Ij 

and represents a transformation y » ax+b of the real line. 

Mow we can discuss the question of obtaining solutions to the 
equation 


( 8 ) 


{ i - ( i-o “ii-i ) i. ' 


j.(0) = x Q e k" , 

It is well known (17], section 3) that the solution to (8) is given by 

(9) I x(t) = $(t)x^ 

where $(t) is an nxn matrix, called the transition matrix for (8) , which 
is the solution of 





I 


1 


1 
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£(t) “ (aq + I J i , *(0) - I . 

We are» therefore, interested in solutions of equation (10). Kow it can 
be shown ([8], Theorem 5) that for all t, ^(t) belongs to the matrix Lie 
Group {exp{A .} that is, the matrix 4>(t) will evolve on the Lie Group 

associated with the Lie Algebra generated by the coefficient matrices 



There are two important consequences of this fact for power conver- 
sion networks described by bilinear equations. The first of these stems 
from the fact that in addition to having the properties of a group, a 
Lie Group has the properties of a manifold, that is, a subset of Euclidean 
space IR^ with special geometrical characteristics. The two-dimensional 
surface of a sphere in three-dimensional space is an example of a manifold. 
The geometrical characteristics of the particular manifold on which the 
state of a bilinear system evolves will play a fundamental role in deter- 
mining the nature of the behavior of the system. For example, if the Lie 
group is compact, i.e. closed and bounded as a subset of then we know 

that the state of the system is bounded, i.e. the amount of energy stored 
in the inductors and capacitors will be finite. For instance, the Lie 
group S0(n) is bounded, and the group Sp(n) is not. Thus, given a bilinear 
system of equations, the natural first question we ask is "On what Lie 
group does its state transition matrix evolve?" 

The second consequence is that questions about the system (such as 
controllability, reachability, observability, and stabilizability) can be 
reduced to questions about the Lie algebra generated by the coefficient 
matrices. In many cases it is possible to arrive at conditions which are 
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easily visualized and tested, as for example our results of Chapter 4, or 

the results of reference (6] . Of relevance here is the classification 

result from the theory of Lie algebras, by which any simple Lie algebra is 

shown to be equivalent to one of a short list of canonical algebras, 

([28] Chapter 2, [24]). For solvable algebras there is no such complete 

list of all possibilities, although a good many facts are known, [28] . 

Note that if is a nonsingular matrix, equation (8) is unchanged by the 

change of variable x "*■ 2L equation (10) is unchanged by 

Now let us further discuss the problem of obtaining a solution to 

equation (10). Let , Xj ® basis for the Lie algebra L 

generated by {A , A ,*««,A } . In reference [19] Magnus showed that 

•“0 1 — in 

there exists a t^ such that for 0 $ t < t^ the solution to (10) can be 
expressed in the form 


( 11 ) 


{ i 


(t) 


b^(t)X^+ b^CUXj +---+bj^(t)4 


where bj^,***,bj^ are scalar functions of time satisfying differential 

equations which depend on L and u, ,***,u . The difficulties with this 

i. m 

approach to obtaining a solution are that it is difficult to derive and 
solve the differential equations for bj^,***,bj^, and that only in severely 
restricted cases [31] is the representation (11) valid for all t in [0,®). 

Wei and Norman [31] showed that it is often preferable to look for a 
solution in the form of a product of exponentials. In fact they showed 
that there exists a t^ such that for 0 4 t < t^ the solution to (10) 
can be expressed in the form 


I i . 
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U2) I *(t) 


g^uix^ 

e e • • • e 


where scalar functions of time satisfying differential 

equations which depend on L and on u_,***u • Moreover, if L is solvable, 

1 m 

or if the matrices are 2x2, then this representation is global, i.e. 
it holds for all 0 4 t < . In § 3.5(c) we give an example of such a 

global representation. 

In conclusion we note that while in some cases it may be possible to 
analyze the behavior of electrical power conversion networks with piecewise 
linear models, concepts from the theory of Lie groups and Lie algebras are 
useful in characterizing the inherent dynamical features of such systems, 
especially since the methods and conclusions are basis-free, i.e. they 
do not depend on the particular basis chosen for the state space. We 
give examples in § 3.5. 


§ 3.5 Network Examples 
(a) An SQ(3) Network 

Pig. 3.8 shows a simple network in which charge stored on one capacitor 

can be transferred by means of the inductor to the other capacitor. If 

the capacitances are different a voltage conversion will be effected. A 

transferral cycle might be executed as follow’S. Starting from 

V (0) a V , 1 (0) a 0, V-(0) e 0, thc switch is held in the u « 0 position 
1 10 3 ^ 

until - 0 and all the energy is stored in L^, at which time the switch 
is changed to the u = 1 position. It is held there until again becomes 
zero, at which time the switch is reverted to the u » 0 position. All the 
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C, i V, 


C2 'T' Vg 


Fig. 3.8 

C, 

initial energy is now stored in C«, with V_ » - r** V . While this network 

• • V2 xo 

is too simple to model a complete DC-DC convertor, it may very well model 
the conversion portion of a convertor in which the charge on is 
replenished from an external supply during the u « 1 portion of the cycle, 
and a load current is drawn from during the u « 0 portion. The egua-- 
tlons for this network are 


CiV^ - 

(l-u)!^ 


- 

“ ^3 


• 

LI » 
3 3 

-(l-u)Vj^ - uVj 

• 


X X 

Letting V. ■ — ^ for i ■ 1, 2 and I ■* — 2. so that •|(x'x) is the total 
stored energy, we obtain 


0 0 
0 0 


“'1 P 

0 + u 1 0 


0 Hi), 


j-W. 0 0 

L ^ 
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where co^ 


•^1 


w. 


1 




3^2 




X ■ |^(1-U)A^ + ^ 2 ]* 


where 

^1“ 

0 . 

0 

“^l 

^2“ 

— 

0 

0 

0 



0 

0 

0 


0 

0 

(0 

2 



-to 

|_-u»i 

0 

0 


0 

-U) 

2 

0 


We see that and A^ are just scalar multiples of and ^ of the Lie 

algebra Example (5) in § 3.4. Thus the Lie algebra generated by A^ and 

A^ is o(3), and the transition matrix £^(t) for this network evolves on the 

2 

Lie group S0(3) . The state vector >c(t) evolves on the 2-sphere S . 

(b) Simple atnd solvable Parts 

If we add a current sink in parallel with of Fig. 3. 8, to represent 
a load for instance, we obtain the network of Fig. 3.9. 

U 



Fig. 3.9 
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Th« State equations are now 



Letting ■ 1 we can use the method of I 3.2 to put this in bilinear 
form« thus: 



i.e. X ■ (B^ uBj^lx 


The Lie algebra L generated by and B^ is six-dimensional and has as a 
basis the following matrices t 


0 

0 



0 0 0 

0-10 
10 0 
0 0 0 


0 

0 



0 10 
0 0 0 
0 0 0 
0 0 0 


0 

1 




-1 

0 

0 

0 


0 

0 

0 

0 



0 








We see that ® reducible representation of o(3) , while 

X^, represent a solvable Lie algebra. ?h-.!s L is decomposable into 

simple and solvable parts, with the simple part determined by the "natural 
dynamics" of the network, i.e. the interconnection of its inductors and 
capacitors, while the solvable pare is contributed by the "driving forces", 
i.e. the batteries and current sources. This type of deccmiposition is 
a general characteristic of the types of networks we are considering, 

(as introduced at the beginning of § 3.2) . 

(c) A Transformerless DC-DC Convertor 

If we now replace the capacitor Cj^ of Fig. 3.9 by a battery and allow 
the switch to take a third position in which the inductor is unconnected, 
we obtain a model of a simple DC-DC convertor, as shown in Fig. 3.10. 



Fig. 3.10 
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He have used two switches in Fig. 3.10 to represent the switching action; 

1 denotes that switch i is closed « and ■ 0 denotes that switch i 
is open* where i « 1, 2. To ensure that the model is well-posed (as 
discussed in § 3.2) we require (a) that u^(t)U 2 (t) - 0 for all t« (i.e. 
both switches cannot be closed simultaneously ) , and (b) that if one switch 
is open, the other c&nnot be opened unless = 0. The conversion cycle 
we envisage is similar to that described in example (a) above, i.e. first 
we set u^ 1 until reaches some predetermined desired value, then we 
let u^ *■ 1 until is again zero. In order to obtain a smoother output 
voltage it may be desirable to use a low-pass filter at the output, as 
depicted in Fig. 3.11. 



Fig. 3.11 

The network of Fig. 3.10 might be implemented with a scheme such as that 
shown in Fig. 3.12, in which ® 1 when transistor T^^ is turned on, and 
U 2 ® 1 when the current in diode is nonzero. 

We shall now illustrate the method of Wei and Norman for obtaining 
solutions to the state- evolution equations of the network of Fig. 3.10. 
These are: 


T 


T 


T 



I 


Fig. 3.12 


0 0 


0 0 B 


-y + u^O 0 O + u^a 0 


0 0 


0 0 0 


0 -a 




0 > X 


0 0 0 I X 


.here Xj - . x^ - = 1, 0 - ^ , 6 - , Y " ^ . 

i.e. ^ X = + U 2 A 2 )x , 

/vft C 

We find that A^] = 0, [A^^, A^] = — — [A^i ?]^] = “^—1 • 

Thus the Lie algebra generated by A^^. A^ is solvable and has A^, 

A as a basis. By the Wei-Norman result we therefore know that there 
-2 

exirrt functions g^Ct) , 9 j^(t) , ^ 2 ^^^ ^ ' 

9tA) ’A 

$(t) = e e e 

where *= + ^1—1 ^ “ i. • 
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To obtain the differential equations satisfied by g^ vre look at 


g g,A, g«A_ g^A^ g,A, g.A_ g«A^ g,A, g,«A_ 

* " 9q^® e e + e g A e e + e e 92^2® 


9q2^ * 9j^ ^ e Aj^e / i + 02 \ ® ® -2 ' - * 


We now make use of the Baker-Hausdorff leuuua (§ 3.4) to obtain: 


” 0 - 0 . 0-0 


. -^Aj. — - Aj . 9oIAj,,Ajl lio'fio'il” * 


»1 • 


Similarly 


’A. 

e A^e 


9 / a ^ “g«A_ 

0 - 0 . ” 0-0 
e A^e 


^•(f) 

‘.-fr) 


We therefore obtain 


i ■ [(^O *f-V2)io * (% -^Vz)^! *"^2 i 

which, on comparison with the defining equation for yields 
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. ^ ctB . 

^0 Y 


« 1 


- h: 


r V2 


= u. 


Since ^(0) = ^ we have 9 q( 0) = =92^0^ * 0. We can therefore write 

the defining differential equations for g^, g^^, g^ as 




r 

1 

0 

1 

Y ^1 


1 


'vo.' 


0 

• 

^1 

= 

0 

1 

n 

B 




9j^(0) 

= 

0 

^2 


0 

0 

1 

> 


'^2 

t 

g2(0) 
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(d) Two Fourth-Order Lossless Networks 

Lossless electrical networks are of interest to us since they may 
define that part of a power conversion network corresponding to the simple 
part of the Lie algebra, as discussed in example (b) above. Here we show 
how a small change in network topology yields a fundamental change in tho 
associated Lie algebra. Consider the network of Fig. 3.13, where we assume 
that the two switches are operated synchronously, i.e. the single control 
variable u denotes the state of both switches. As usual we let = ^1*^' 
X- = V_*^, etc. , and we let a * , B * — ^ , Y * — — ■ » 6 = — ^ . 

■^2 ’SS -^4 '^2 

Then we obtain 


— t 


1 


-.4 . 
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l^s 
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Fig. 3.13 


where 


-0 




0 -a 

0 o' 


0 0 0 -Y 

a 0 

0 0 

K = 

0 C 5 0 

0 0 

0 -B 

-1 

0 

1 

o» 

o 

o 

0 0 

B 0 


o 

o 

o 

J 


Now the Lie algebra generated by ^ and is six-dimensional, except 
when o “ B or Y “ 6 when it is four-dimensional. If a « B and y = 6 then 
a B B " Y * ^ the Lie algebra is two-dimensional, with basis A^^. 

NOW the network of Fig. 3.13 can be redrawn as shown in Fig. 3.14. 



Fig. 3.14 


We can modify this network slightly by adding another synchronized switch 
which has the effect of reversing the polarity of C^. Fig. 3.15 depicts 
this siti'.ation. 
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1 


1 



with a, 6f Y* ^ defined as before. Again we find that the Lie algebra 
generated by ^ and is six-dimensional, except when a “ B or y = 5 
when it is four-dimensional. This time however, ifa“6“Y®5 then the 
Lie algebra is three-dimensional, and in fact is a representation for the 
Lie algebra o(3) . Thus, when a = 0 « y = ^ the state transition matrix for 
the network of Fig. 3.15 evolves on the Lie group S0(3) . 


Higher-Order SO (3) .Networks 


The abstract Lie algebra o(3) is defined by the relationships 


S , [S , S ] - S 
-X -a “X -y 


IS , s ] = s 
-X -y -« 


A representation of this abstract Lie algebra is a set of three matrices 
which satisfy these relationships. A representation is said to be 
irreducible if its component matrices cannot simultaneously be put in block 
triangular form, as in § 3.3. Now it can be shown ([28], Chapter 1) that, 
over the complex field, all irreducible iixn irreducible representations 
of 0(3) are equivalent to the following representation, where J = and 
® i(n-i) : 
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^-1 « 


%-T 


0 0 
0 

0 -»ir 0 


° ’^-1 

"^n-l ° 


S - /^ 


An nxn complex representation of o(3) can be made into a 2n x 2n real 
representation by identifying a given complex matrix (R + 2 ) with the 

real matrix 
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£ 1 


For instance* when n « 3 we obtain the 6x6 real representation of o(3) 
given by 


A 


— X 


1 


• ^ 



1 0 -1 0 


0 1 0 1 

1 

1-1 0-1 

A - ± 

-1 0 1 I 

1 

/2 

1 

1 

1 0 -1 0 


0-101 



0 1 0 j 


1 

1 0 1 0 
1 

< 

1 0 11 
1 


|-1 0 1 

1 

0 1 0 I 


1 

, 0-1 0 






Suppose now that we would like to find a sixth-order network which has the 
state equations 


X. [> 


uA + (l-u)A 
“X ~y 
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To do this we exendne A atnd A and reverse the process by which we 

-X -y 

normally obtain the state equations for a network. Assume that is an 
inductor current. We see from A that when u ■ 1 we need element 5 to be 

—X 

a capacitor C_ and element 3 to be an Inductor L., with L, and L. connected 
o 3 X J 

in parallel with C_. From A we see that when u = 0 we need L, and L. in 
5 -y 13 

parallel with a capacitor C 2 « And in parallel with inductors and L^. 
The resulting network obtained by this process is shown in Fig. 3.16 in 
which the switches all operate synchronously: the switches denoted — — 
are closed when u ■ 1 and open when u ■ 0, while those denoted -—9^*— 

are closed when u >> 0. and open when u « 1. 

By a similar process we can construct a network of order 2n for 

any n ^ 2 whose state transition matrix evolves on the Lie group SO (3) . 


(f) A Reducible Network 


Our final example shows how reducibility of the Lie algebra can 
correspond to reducibility of the network. In Fig. 3.17 we wish to trans- 
fer energy from the battery E to the output capacitor C.. Me might ask 
whether an unlimited amount of energy can be extracted from the battery. 
The two switches are operated synchronously. In a similar manner as 
before we obtain the state evolution equations as 


X « + uAj^) x 




'•'T 



1 



The Lie algebra generated by and is o(3) , and thus we conclude that 
this network is dynaaiically similar to that of Fig. 3.8 and that only a 
finite amount of snergy can be stored in C^, since the associated 
Lie group is bounded. This may at first seem surprising, since 
the network of Fig. 3.17 has a battery in it. In the case *= » 1 

the Lie algebra has as a basis 
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where - X^^, IX^.X^J - X^, [X^^X^J - X^ . 

Examination of the network of Fig. 3.17 shows that we can reduce 
it to a simpler foirm. The network is unchanged if the switch u and capaci- 
tor are interchanged, giving the network of Fig. 3.8 but for a battery 

in series with C^. Now the total energy which can be extracted from the 

1 2 

combination of Fig. 3. 18 is— CE . (This is an interesting singular 



Fig. 3.18 

optimal control problem which may be solved by observing that the arrange- 
ment of Fig. 3.18 is externally identical with that of Fig. 3.19, since 
both are governed by 

V(t) « E + -^ /* I(T)dT . ) 

*'o 

In fact we see that the network of Fig. 3.17 is essentially the network 


of Fig. 3.8, but with the capacitor being given an extra initial voltage 
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■o I(t) 


V(0)*E 


Fig. 3.19 

of E. The matrices X^ given here are a reducible representation 

of the Lie algebra o(3) : they can be reduced by eliminatirtg from each one 
the last row and coluiw> which represents the contribution of the 
batteries and current sources in the network ^ as wo saw in example 3.5(b). 
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CHAPTER 4 

FEEDBACK STABILIZATION OP BILINEAR SYSTEMS 


In Chapter 3 we showed that bilinear systems arise naturally in the 
study of power conversion networks. Since it will in many cases be desired 
to stabilize such systems by means of state feedback, we are interested in 
the general study of feedback stabilization for bilinear systems. In 
Chapter 1 we discussed this question for a certain class of DC-DC conver- 
sion networks. In Chanter 4 we shall address ourselves to the question of 
feedback stedsilization of systems which evjly/e on the (n-1) sphere S*' 
i.e. systems whose state vector x e IR*^ satisfies x’x = constant. We have 
in mind two specific questions concerning the feedback stabilization of a 
general bilinear system: (i) If the system is controllable (in some 

approximately defired sense ) , can we find a feedback law such that the 
closed-loop system is asymptotically stable about a particular point? 

(ii) If the system is controlleible, can we find a feedback law such that a 
given oscillation is stabilized for the closed-loop system? 

Theorem 4.1 Consider the system of equations 

m 

k~ 1 ^ » X* (0)£ x(0) = 1, 

i=l 

where £A^+ Al£ = £ for 1 i $ ro and £' = Q > £. Given some 3^ which 

satisfies suppose that the matrix whose columns are ' 

QA.x.,***, QA X- has remk n-1. Let u = (u, u_ • • • u ] . Then there 
A. —2 — U ^ “tU “ U — 1 7 m 

exists a feedback control law u = 2i(x) such that the closed- loop feedback 
system is asymptotically stable from any starting point x(0) other than -Xq» 
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Proof. Let M denote tl»e set of points x satisfying x* J^x « 1; when 2. “ i 

then M is the (n-1) sphere s” First, note that we have x^' (t)^ ^(t) =1 for 
all t in [0, ») if x(0) e M, since ^ .^.'2.2. “ i'22 2*22 ® ® since 

22i + A^2 » 0 for each i. Next, consider the Lyapunov function 

V = j ^2“2o^' 2 ^2“2 q^‘ Since g > 0 we know from the Lyapunov theorem 

of §1.4(e) that if in some subset N of M with have V $ 0 with 

• 

V a 0 only along the trajectory x(t) 2 q , then the desired asymptotic 
stability about x^ will be obtained. 

Now since x(t) £ M for all t we have 

V = ^ x'2(2 “ 2 q^ ■ 2.0^ ' £2 

thus the desired stability about XqS Ncm will ensue if 21 'q 22 5 0 in N, 
with “ 0 only at x^ . 

Now let f(o) be any odd function on the real line for which f(C) =0 
only at * 0. Let u^ = f(x,*Q22i2^ 1 '■ i < m. Then 

2’o 22 = 2*0 2 I f(2’oa2i2)2i2 

is 1 

m 

“ I (2'n22<2>^<2'nS2H2J 

i^l u 1 U X 

m 

= I (2'A:22o’f<s'A:22o) 

i=i 

m 

“ I ^?l'22i2o ^ ^ ^2'22i2o^ since £A^ + = 0 and f is 

odd. 


1 


J 


i 





1 


t 


-165- 

Thus ^ ^ *11 X € M, and x.'q2^x = 0 if and only if x’g^ ^ 

for 1 $ i $ n. Now since tQ A ,x ^ ,0 A _x . ,♦ • • ,0 A x ] is of rank n-1 and 
Jt^is in the one- dimensional subspace perr^endicular to ^2^]^2 Lq ' * ** '2.— m— 0 ^ ' 
we must have that 

x'QA.x^ = 0 for all i x = kx . 

for some real k. From x'2.x = 1 we obtain k = 1, and the result follows. 

QED 

The following theorem was proved jointly by Professor R. W . Brockett and 
myself. 

Theorem 4.2 . Consider the system of equations 

X = (A‘+ uB)x , x(t) eiR^, x' (O)jt(O) = 1, 

where A' + A = B* + B = 0^. Suppose that Xq^Iq “ 1 that A = 0 and 
B^jc Q / 0. Let u = where f (a) is any function on the real line 

satisfying of (a) ^ 0, with f(a) =0 only at a = 0. Then the resulting 

closed- loop system is asymptotically stable about in a neigliborhood N 

of Xq if and only if the pair ( 21 'q§»A) is observable (in the linear system 
sense) . Moreover, a sufficient but not necessary condition for this is that 
the set of matrices {Ad B^a} span the space of skewsymmetric matrixes, 
where k » and Ad is as defined in §2,4. 

Proof . Since x = Ax + f(x'QBx)B x we have x'qX = (x'qB x) f ^ ^ 

for all X. Now x'^x = 0 for all t if and only if x'Q®ii “ ® 

At 

i.e. if and only if x'q® ® t for some x^« Thus, by the 

Lyapunov method in the proof of Theorem 4.1, we conclude that asymptotic 

stability in N is equivalent to the requirement that there does not exist 

At 

an Xj^ ^ Xq in N such that x'qB ® “ 2ii ® ° t.But this is equivalent 
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'i 
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( [7 ] , section 13) to the requirement that there is no in N such that 

“ 0 for i ■ which is equivalent to the requirement 

that be observcO^le (17], section 14). Thus, asymptotic stability 

ensues if and only if (x'qB,A) is observable. Now if x'^ A * 0, it is 

straightforward to show by inauction that x' Ad!”B= (-1)"* x* BA™ for 

—0 A— —0—— 

(E*g: (A B - 2. A) = - x*q 2. ^ ^ * Thus asymptotic 

stability in N ensues if there is no x. in N such that x* Ad.^B x, - 0 

—X o A 1 

for all k and x'q A X| » 0. But x'q Ad^^^B x^^ = tr Ad^^B , anc x'^ 

•= tr A x^x'q. Now if tr x y* * 0 for all X satisfying X + i’ “ 2. ' then 

Y = y * } and if jc is symmetric then ^ ® c x for some c e IR. Hence, if the 

matrices (Ad^^B^.A) speui the skewsymmetric matrices then x'^Ad^^B x^^ = 

i'o — — 1 “ ® will iit^)ly that x^^ “ i 21 q (assuming that il^ » 21 q ^ ’ 

i.e. If {Ad^^^,A} spans the space of skewsymmetric matrices then asymptotic 

n"*l 

stability ensues, from any starting point in S other than To show 

that thiis condition is not necessary, consider the case where 


0 10 0 0 
-10 10 0 
0-1010 
0 0-100 
0 0 0 0 0 


0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 1 

000-10 


Then the Lie algebra generated by A and B^ is the space of 5 x 5 skewsymmetric 
matrices, but {Ad^ matrices of the form 


0 1 0 0 d 

-10 10c 

0 -1 0 1 b 

0 0 -1 0 a 

-d -c -b -a 0 
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whlch Is not the ^ole space of 5 x 5 skevrsymmctric matrices. Now the only 
in S** ^ satisfying A x^ « £ is 



Prtxn this we find that B, A) is observable, and the result follows. 

QED 

In general, it is not possible to obtain a global stabilization on 
^ for systems of the type considered in Theorems 4.1 and 4.2 when the 
feedback control law \j(x) is restricted to being a continuous function. 
(Cf. Hopf's Theorem concerning the number of singular points of a smooth 
vector field on a manifold without boundary) . In practice the fact that 
there will always be a "deadpoint" x such that x|^ » 0^ would pro- 
bably be of little concern. However, one might ask the question "Can the 
point Jt^ be chosen to be any other point on S ?" We would expect that 
it can, since this amounts to a smooth topological deformation of the 
vector field obtained in Theorems 4.1 and 4.2. As an example, consider 
the system on given by 



Suppose that it is desired to stabilize this about the point (a,b) where 
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2 2 

a +b =1 and a,b > 0. Then let u = a •• Xj^. One can show that asymptotic 
stability about (a,b) ensues, from any starting point in other than 
(a,-b) . The vector field is as shown in Figure 4.1. 



Fig. 4.1 


A similar example can be given on S . In such cases one proves stability 
in s“ -{d} about £ by showing that there is a neighborhood about £ 
with £*£ > 0 in N^, and that there is a neighborhood about £ with 

d'x ^ 0 in N^, such that U * s" ^ - {d}. 

Finally, we consider the problem of stabilizing an oscillation on 
S*' Given £ = f (x.,£) one should choose, if possible, £(j0 so that 
X E « {x|v(x) = o} where V is a suitable Lyapunov function and E is the 
set of points in the desired orbit. One should also choose u(x) so that 
on E £ follows the orbit cyclically. 

2 

He consider first an example on S . Suppose 
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and that it is desired to stabilise an oscillation around the set {x^ s a } , 
for some 0 ^ a < 1. 







Fig. 4.2 


Let «2 * that when Uj^ * 0 on {xj^ 


a} we have the simple harmonic 


oscillation given by 


: ':ll: 


It then remains to choose u^(x) so that x^ a. 

We shall make use of the Invariance Principle of LaSalle, [17 ] , 


[371. Consider the periodic or time-invariant system of equations 
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! x(t) ■ f,(x,t) • f(x,t+T) 
x(t) e r” ; 0 <t«»} jc(0) " Xq • 


Let V(3c,t) e V(x,t+T) be a Lyapunov function (not necessarily positive 
definite) on a closed bounded set G for this system of equations, as 
defined in § 1.4(e). Let E = {^|w(x^) * 0, x e g} where V(x, t) W(x^) $ 0. 
Then the Invariance Principle of LaSalle (from which our Theorem 1.3 
follows) states that any solution of the above equations which remains in 
G for all t ^ 0 approaches a subset of E which is the union of all the 
invariant sets which lie entirely within E. (An invariant set H is one 
for which x,(tj^) e H ^ ^2 ^ ^1^ * 

Now since x, - let us try Uj^ = Then x^^ » (a-Xj^)x 

2 2 

Consider the closed bounded subset N^ of S defined by N^={}c e S | x^ 4 a} . 
Let where a = (a, 0, 0). 

Then » “*.'x 

/ X 2 

- -a(a-Xj)x2 
0 in Nj^. 

Now Xj^ * a defines a trajectory for the system, thus no other trajectory 
can cross {x^^ = a} , i.e. Any motion starting in N^^ at t=0 remains in N^^ 
for all t » 0. 

Let « {x e Nj^jVj^(x) * 0} 

2 

* e S I Xj^ i a and either x^ “ 0 or Xj^ ■ a} . 
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On {x^ “ 0} we have ■ “X^ ^ 0 on Ej^ only at (0, 0, -IJ . Thus the union 
of the invariant sets in is 

{x e Nj^jx^ * a or 3^' ■ ( 0 , 0 , - 1 )}. 

Furthermore a motion starting from any point in other than (0, 0, -I) 
cannot approach (0, 0, -1), since x^^ ^ 0. Thus, by the Invariance 

Principle, we conclude that the desired oscillation is stabilized from any 
starting point in other than (0, 0, -1) . 

2 2 

Considering the subset of S defined by * (x e S |x^^ ^ a} 
together with •* a'x , we conclude by a similar argument that the oscill- 
ation is stabilized from any starting point in other than (0, 0, 1). 

Thus, we have stabilized a circular (simple harmonic) oscillation on 
2 

the sphere S around Xj^ “ a from any starting point other than (0, 0, ±1) . 

Furthermore this means that we have stabilized such an oscillation around 

2 

any circular orbit on S , since bilinear evolution equations are invariant 
under the transformation x x P , as mentioned in § 3.4. 




i- 


« 


j 



! 




i 


\ 


1 


I 


REFERENCES 

(1] B.D.O. Anderson and J.B. Moore, Linear Optimal Control , Prentice-Hall, 

New Jersey, 1971. 

[2] N. Balabanian, Network Synthesis , Prentice-Hall, New Jersey, 1958. 

t31 J.G. Belinfante, B. Kolman, and H.A. Smith, "An Introduction to Lie 

Groups and Lie Algebras with Applications", SIAM Review ; Vol. 8, Ho. 1, 
Jan. 1966 (Part I); Vol. 10, No. 2, April 1968 (Part II); Vol. 11, 

No. 4, Oct. I960 (Part III). 

(4) R. Bellman, Stability Tlieor y o f Piff erential rn u ations , McGrav; Hill, 

New York, 195liri)ovir,"New York) 1969. 

(5) R.W. Brockett, "Poles, Zeros, and Feedback: State Space Interpretation", 
IEEE Trans, on Automatic Control , Vol. AC- 10, No. 2, April 1965, 

pp. 129-135. 

16] R.W. Brockett and J.L. Willems, "Frequency Domain Stability Criteria: 

Part I", IEEE Trans. Automatic Control , Vol. AC- 10, pp. 255-P61, 1565. 

(7] R.W. Brockett, Finite Dimensional Linear Systems , Wiley, New York, 1970. 

[8] R.W. Brockett, "System Theory on Group Manifolds and Coset Spaces", 

SIAM J. on Control , Vol. 10, No. 2, May 1972, pp. 265-284. 

19] R.W. Brockett, "On the Algebraic Structure of Bilinear Systems", in 
‘Theory and Aop l icatiens of Vari>ii,lG Structure Systems , R. Mohler 
and A. Ruberti Eds., Academic Press, New York, 1972. 

(10] C. Chevalley, Theory of Lie Groups , Princeton University Press, 

Princeton, 1946. 

(11] P.M. DeRusso, R.J. Roy, C.M. Close, State Variables for Engineers , 

Wiley, New York, 1965. 

(12] P.L. Falb and G. Zames, "Or. Cross-Correlation Bounds and the Positivity 
of Certain Nonlinear Operators" and "On the Stability of Systems with 
Monotone and Odd-Monotone Nonlinearities", IEEE Trans, on Automatic 
Control , Vol. AC-12, pp. 219-223, April 1967. 

(13] A.F. Filippov, "Differential Equations with a Discontinuous Pight-Hand 
Side", The American Mathematical Society Translations , Series 2, Vol. 42, 
pp. 199-230. 

(14] E.R. Hnatek, Design of Solid State Power Supplies , Van Nostrand Reinhold, 
Nev? York, 1971. 

(15] F.F. Judd and C.T. Chen, "Analysis and Optimal Design of Self-Oscillating 
DC-to-DC Converters", IEEE Trans, on Circuit Theory, Vol. CT-18, Mo. 6, 
Nov. 1971, pp. 651-658. 


PRECEDING PAGE BLANK NOT FILME2> 


- 173 - 







1 


I 


[16] S. Kami, network Theory; Analysis and Synthesis, Allyn and Bacon, 

Boston, 1966. 

[17] J.p. LaSalle, "An Invariance Principle in the Theory of Stability", in 
Differential Equations £uid Dynamical Systems , J.K. Hale and J.P. 

LaSalle Eds., Academic Press, New York, 1967, pp. 277-286. 

[18] S. Lefschetz, Stability of Nonlinear Control Systems , Academic Press, 

New York, 1965. 

[19] W. Magnus, "On the Exponential Solution of Differential Equations for 
a. Linear Operator", Comm. Pure Annl. Math, Vol.7, 1954, pp. 649-673. 

[20] L. Markus, "Dynamical Systems on Group Manifolds", in Differential 
Equations and ! >/nanical Systems , J.K. Hale and J.P. LaSalle Eds., 

Academic Press, New York, 1967, pp. 487-493. 

[21] J. Meixner, "On the Theory of Linear Passive Systems", Archive f or 
Rational Mechanics and Analysis , Volume 17, 1964, pp. 27R-296. 

[22] A.N. Michel and D.w. Porter, "Practical Stability and Finite-Time 
Stability of Discontinuous Systems", IEEE Trans, on Circuit Theory , 

Vol. CT-19, No. 2, March 1972. 

[23] E.T. Moore and T.G. Wilson, "Basic Considerations for DC to DC conversion 
networks", lEEi: Trans, on Magnetics , Vol. 2, No. 3, Sept. 1966, 

pp. 620-624. 

[24] S. Murakami, "Sur la Classification des AlgSbres de Lie Reelles et 
Simples", Osaka J. Math . , Vol. 2, 19C5, pp. 291-307. 

[25] J.R. Nowicki, Power Supplies for Electronic Equipment , Volumes I and II, 
Leonard Hill, London, 1971. 

[26] R.P. O'Shea, "A Combined Frequenc'y-Time Domain Stability Criterion for 
Autonomous Continuous Systems", ILEE Trans, on Automatic Control, V ol. 
AC-li, pp. 477-483, July 1966; Proc. JACC, pp. 832-840, 1966. 

[27] W. Rudin, Principles q 6 Mathematical Analysis, McGraw-Hill, New York, 
1964. 

[28] H. Samelson, Notes on Lie Algebras , Van Nostrand, New York, 1969. 

[29] E.C. Titchm.arsh, Introduction to the Theory of Fourier Integrals , 

Oxford University Press, Oxford, 1937. 

[30] F.W. Warner, Foundations of Differentiable Manifolds and Lie Group s, 

Scott Foresman, Glenview, Illinois, 1971, 

[31] J. Wei and E. Norman, "On Global Representations of the Solutions of 
Linear Differential Equations as a Product of Exponentials", Proc. 

Am. Math. Soc. , April 1964, pp. 327-334. 


i 


J 


i 



- 175 - 


(32] J.C. Willems and M. Gruber# "Comments on "A Combined Frequency-Time 
Domain Stability Criterion for Autonomous Continuous Systems"", 

IEEE Trans, on Automatic Control, Vol. AC-12, pp. 217-219, April 1967. 

(331 J.C. Willems and R.W. Brockett, "Some New Rearrangement Inequalities 
Having Application in Stability Analysis", IEEE Trans, on Automatic 
Control , Vol. AC-13, No. 5, Oct. 1968, pp. 539-549. 

(34) J.C. Willems, The Analysis of Feedback Systems , M.I.T. Press, Cambridge, 
Mass., 1971. 

(35’ J.C. 'Willems, "Dissipative Dynamical Systems, Pari; I: Ccncral Tiiccry", 
Archive for Rational Mechanics and Analysis, Volume 45, No. 5, 1972, 
pp. 321-351. 

(36] J.C. Willems, "Dissipative Dynamical Systems, Part II: Linear Systems 
with Quadratic Supply Rates", Archive for Rational Mechanics and 
Analysis , Volume 45, No. 5, 1972, pp. 352-393. 

(37] J.L. Willems, Stability Theory of Dynamical Systems , Nelson, London, 
1970. 

(38] D.H. Wolaver, Fundamental Study of DC to DC Conversion Systems , Ph.D. 
Thesis, M.I.T. , Cambridge, Mass., 1969. 

(39] D.H. Wolaver, "Basic Constraints from Graph Theory' for DC-to-DC Con- 
version Networks", IEEE Trans. Circuit Theory , Vol. CT-19, No. 6, 

Nov. 1972, pp. 640-648. 

(40] T. Yoshizawa, "Asymptotic Behaviour of Solutions of a System of 
Differential Equations", Contributions to Difforontial Equations , 

Vol. 1, No. 3, 1963, pp37i-387. 

(41] D.C. Youla, L.J. Castriota, H.J. Carlin, "Bounded Real Scattering 
Matrices and the Foundations of Linear Passive Network Theory", 

Trans. IP£ Circuit Theory CT-4, 1959, pp. 102-124. 

(42] G. Zames, "On the Input-Output Stability of Nonlinear Feedback Systems, 
Parts I and 11", IEEE Trans, on Automatic Control , Vol. AC-11, pp. 
228-238 and 465-476, April and July 1966. 


I 



